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Abstract 


As calculations of many-body nuclear systems become more and more precise, 
the question has shifted from “can we calculate observables?” to “how precise can 
onr calculations be?”. Improvements in many-body methods have identified the need 
for better few-body potentials as input. There are a wide variety of input poten¬ 
tials, which all accurately describe few-body systems, and all must be transformed 
with similarity renormalization group (SRG) transformations before they are com¬ 
putationally efficient enough to be used in certain many-body calculations. It was 
realized that modern realistic 2-body potentials with very different matrix elements 
evolve under SRG transformations to the same universal low-energy shape. Under¬ 
standing the requirements for universality in evolved potential matrix elements can 
aid in the construction of better potentials for the many-body problem. Furthermore, 
understanding the precision of few-body evolution is a key step to setting error bars 
on theoretical predictions. 

We first examine how the universality of two-nucleon interactions evolved us¬ 
ing similarity renormalization group (SRG) transformations correlates with T-matrix 
equivalence, with the ultimate goal of gaining insight into universality for three- 
nucleon forces. Because potentials are fit to low-energy data, they are (approxi¬ 
mately) phase-shift equivalent only up to a certain energy, and we find universality in 



evolved potentials up to the corresponding momentum. More generally, we find uni¬ 
versality in local energy regions, reflecting a local decoupling by the SRG. The further 
requirements for universality in evolved potential matrix elements are explored using 
two simple alternative potentials. We see evidence that in addition to predicting the 
same observables, common long-range behavior (i.e., explicit pion physics) is required 
for universality. In agreement with observations made previously for V\ ow k evolution, 
regions of universal potential matrix elements are restricted to where half-on-shell 
T-matrix equivalence holds. 

To continue the study in the 3-body sector, we create a simple 1-D spinless bo¬ 
son “theoretical laboratory” for a dramatic improvement in computational efficiency. 
We introduce a basis-transformation, harmonic oscillator (HO) basis, which is used 
for current many-body calculations and discuss the imposed truncations. We con¬ 
firm that evolution to universal low-energy 2-body potential matrix elements is the 
same for 1-D bosons as 3-D fermions, and show that a further simplification of using 
positive-valued eigenvalues rather than phase shifts is valid. When SRG evolving in 
a HO-basis, we show that the evolved matrix elements, once transformed back into 
momentum-representation, differ from those when evolving with momentum repre¬ 
sentation. This is because the generator in each basis is not exactly the same due to 
the truncation. In the 2-body sector, this can be avoided by increasing the basis size, 
but it remains unclear whether this is possible in the 3-body sector, as computational 
power required is greatly increased for three-body evolution. 

In our efforts to study universal matrix elements in the 3-body sector by ob¬ 
serving momentum representation matrix elements, we observe oscillations much like 
those appearing from truncation errors in the 2-body sector. We can identify that 



the spectator particle adds strength in far-off diagonal potential matrix elements of 
the embedded 2-body potential, thus truncation errors appear in 3-body HO-basis 
evolution at much higher values of the decoupling scale than expected from 2-body 
calculations. Observing matrix elements of a part of the SRG evolution that should 
be the same in 2- and 3-body sectors shows that they are different, which indicates the 
difference is an error in the HO-basis 3-body evolution. With a better understanding 
of the source of errors in 3-body input potentials, we hope to gain further insight into 
how to make progress towards precision nuclear many-body calculations. 
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Chapter 1: INTRODUCTION 


The nucleus is made up of positively charged protons and neutral neutrons, and 
thus the force binding it together cannot be electromagnetic in origin. The strong 
interaction that is responsible for this binding has been the focus of study of many 
physicists for over 80 years, since the discovery of the neutron in 1932. Despite many 
triumphs, however, there is still much work to be done to achieve accuracy goals for 
calculations of many-nucleon systems. 

Initially, much of the theory of low-energy nuclear physics was driven by exper¬ 
iment, including the discovery of new particles. Yukawa proposed that the force 
between nucleons is caused by exchanging a particle, later called a pion. This in¬ 
teraction was much like the exchange of massive photons, and with the discovery of 
the pion in 1947, meson-exchange potentials became more sophisticated. In the 50’s 
and 60’s more mesons were discovered and the interaction between nucleons could 
be modeled by exchange of single mesons. In the 70’s more advanced models includ¬ 
ing multiple meson exchange (such as 2-pion exchange) were created. At this time, 
however, it was not possible to perform accurate few- and many-body calculations 
with the potentials due to insufficient computing power, which slowed progress. By 
the end of the 70’s quantum chromodynamics (QCD) was largely established as the 
fundamental theory of the strong interaction, which called into question the use of 


1 



meson-exchange potentials, as mesons are not fundamental particles in QCD. At the 
time, QCD was not solvable in the low-energy regime in which nuclei exist, becoming 
perturbative only at much higher energies due to asymptotic freedom. Because of 
this, meson-exchange potentials were still used, but deemed phenomenological. 

Phenomenological potentials became sophisticated enough to accurately reproduce 
2-body observables and required certain considerations in their construction. Nuclear 
saturation implies that the potential must be short-ranged and strongly attractive 
at separations of a few fm. The deuteron (a nucleus made of one proton and one 
neutron) has an electric quadrupole moment, which means that it is not spherically 
symmetric, and therefore the wavefunction is not purely S-wave, but a mixture of par¬ 
tial waves. This requires a non-central force, which has a tensor character. Scattering 
experiments imply that the interaction turns repulsive at higher energies, which led 
to a “hard core” being put into the interaction [T]. One such phenomenological po¬ 
tential we use in this study is Argonne Ui 8 M- These potentials are very accurate 
in the 2-body sector, but the hard core is prohibitive in the many-body regime for 
most solution methods, because the necessary single-particle basis sizes are too large 
for even today’s computing capabilities. Some of the phenomenological potentials 
lead to triumphs in calculations. First, they can very accurately reproduce 2-nucleon 
observables, and with them irrefutable evidence for the necessity of 3-body forces 
was found. Secondly, they are local in particle separation, which is a requirement for 
quantum Monte Carlo calculations. This is a particular class of many-body methods 
which can accurately describe energy spectra and shapes up to 12 C. 

Around 1990, another breed of potentials began to be developed. Weinberg ap¬ 
plied effective field theory methods to multi-nucleon systems by constructing the most 


2 


general Lagrangian consistent with low-energy QCD The new theory was called 
chiral effective held theory (yEFT), and the degrees of freedom were once again nu¬ 
cleons and pions. yEFT potentials have beneficial qualities; they have the relevant 
degrees of freedom for low-energy nuclear physics, are consistent with QCD symme¬ 
tries, and are much softer than the hard-core potentials. Both phenomenological and 
yEFT potentials are used in modern calculations and continue to be improved upon. 
The most beneficial quality of yEFT potentials is that they are derived in a model- 
independent systematic expansion, which makes theoretical error estimates possible. 
Accurate yEFT potentials are much softer than the hard-core potentials, but they 
still are “too hard” for convergence in large many-nucleon systems. 

To address this problem, Lee-Suzuki transformations were applied in free space to 
integrate out high-energy degrees of freedom and soften an initial realistic potential, 
generating phase-shift-equivalent “low-momentum” or “V[ ow fc” potentials 0E]- This 
can be done in one step or incrementally using a renormalization group (RG) equation 
for the potential |7j. Bogner and collaborators observed that a wide variety of realistic 
potentials have very similar low-momentum matrix elements after softening, which 
they termed the model independence of Vj ow fc potentials Pi- The diagonal V\ ov ,k 
potential matrix elements were found to match in regions of phase-shift equivalence 
of the realistic potentials while the off-diagonal matrix elements matched in regions 
of half-on-shell (HOS) T-matrix equivalence [5]. They suggested that differences in 
the HOS T-matrix and thus the off-diagonal Viowfc potential matrix elements occur 
because of different treatments of pion physics [5|. 

Subsequently, similarity renormalization group (SRG) unitary transformations 
have been used to soften nuclear potentials while preserving observables punumEi 
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H30E|. The Lee-Suzuki transformations are derived using the T-matrix, and exten¬ 
sion into the 3-body sector is difficult, whereas SRG transformations in the 3-body 
sector are conceptually straight-forward. Like Idowfc transformations, the SRG decou¬ 
ples high-energy from low-energy physics, allowing one to truncate the matrices above 
some decoupling scale mmm- Further, the low-energy matrix elements of initial 
realistic potentials also flow to the same form, but differ in detail from those found 
using Viowfc transformations. There is preliminary evidence that the SRG flow to 
common matrix elements extends to three-body forces mm, which are important 
ingredients for consistent treatments of nuclei with RG methods pans]. 

In analogy to the behavior of other Hamiltonians under RG transformations, this 
model independence is naturally interpreted as a flow to universality of the evolved 
potential matrix elements. Operators can be classified according to how their dimen¬ 
sionless coupling changes under RG transformations as relevant (coupling increase 
with decreasing resolution scale), marginal (coupling remains of the same order), and 
irrelevant (coupling decreases) [S|- The universal ability of different potentials to 
describe the same behavior is a consequence of the potentials possessing the same 
relevant and marginal operators, and the flow toward universal potential matrix el¬ 
ements naturally occurs as irrelevant operators diminish. This form of universality 
can have powerful consequences if it can be understood and exploited. It suggests 
that for low-energy problems, a broad class of starting potentials that fit data will be 
nearly equivalent after evolution [201 El]. If realized for many-body forces, it may be 
possible to more easily construct accurate potentials (choosing operators based solely 
on their ease of use, then fitting constants to data), if they flow to a universal form 
after running the SRG. 
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In the effort to improve the nucleon-nucleon interaction, much attention is given 
to improving the potentials from phenomenological considerations or including more 
terms and degrees of freedom in ytEFT. Because the SRG is such an important tool 
to increase computational efficiency, and because the realistic potentials flow toward 
a universal form, in this study we focus on better understanding criteria for univer¬ 
sal SRG evolved potential matrix elements and identifying errors in potential matrix 
elements after few-body harmonic oscillator (HO)-basis SRG evolution. Many of the 
results of previous SRG studies have used 2-body evolution in partial-wave momen¬ 
tum representation [221 ESI El II51H11 El], but in practice, because of the importance 
of the 3-body forces, many-body problems require 3-body potentials and benefit from 
evolution in HO-basis as input (3-body momentum representation evolution has re¬ 
cently been done, but is not widely applied yet pa]). Errors in few-body HO-basis 
evolved Hamiltonians cannot always be distinguished in 3-body observables because 
even the approximate transformations are unitary. These errors can lead to lower 
precision in many-body calculations. Because these errors in the 3-body Hamilto¬ 
nian do not affect 3-body observables but can effect many-body observables, we can 
naturally interpret them as the introduction of spurious 4-body and further forces. 
As unevolved potentials are increasingly more difficult to generate, it is important 
to utilize short-cuts gained from universality in evolved potential matrix elements. 
We also must understand the precision of few-body SRG evolution, which in turn 
will lead to better input Hamiltonians and a better understanding of propagation of 
uncertainties for many-body methods. 
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1.1 Motivation to Improve the Few-Body Hamiltonian 


As mentioned earlier, phenomenological potentials were able to reproduce 2-body 
observables with high accuracy, but attempts at few-body nuclear calculations could 
not reproduce experimental values using these “high-precision” 2-body potentials 
alone. Few-body calculations had enough precision to confirm what many had been 
resisting because of the prohibitive leap in computational difficulty; that interactions 
depending on all three particles are not negligible and must be included. Once de¬ 
veloped, phenomenological 3-body potentials vastly improved the accuracy of bound- 


state calculations for light nuclei. Figure |1.1| shows calculated binding energies of 
light nuclei with only a 2-body potential, and with the inclusion of two different 3- 
body potentials E25I . This picture enforces the idea of a many-body hierarchy; that 
2-body forces are the greatest contribution to observables and 3-body forces are non- 
negligible and must be included for improved accuracy. Higher-body forces can be 
added, each with decreasing contribution. The importance of the 4-body force is still 
an open area of study. 

Modern many-body techniques which use as input modern 2- and 3-body poten¬ 
tials have become very precise. These potentials are able to reproduce light-nuclei 


bound-state energies, but heavier nuclei remain a challenge. Figure 1.2 shows a re¬ 
cent calculation for binding energies and approximated errors of many nuclei using a 
modern many-body calculation [25] . The precision of calculations is a new and de¬ 
veloping held, and to estimate error bands, a calculation is run with several different 
realistic potentials, such that a maximum and minimum value in observables gives 


some idea of the precision. In Fig. L2 we see that within the estimated precision for 
the heavier nuclei, the calculation does not accurately reproduce experimental values. 
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Figure 1.1: Binding energies of various nuclei calculated using a quantum Monte 
Carlo many-body method with only 2-body input potential (AV18), with 2- and 3- 
body input potential (AV18+IL7), and experimental values. Including the 3-body 
shifts the energies by a fraction of the 2-body results, increasing the accuracy system¬ 
atically for all energy levels. Accuracy generally decreases in heavier nuclei. Figure 
from Ref. [25] 
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Figure 1.2: Binding energy per particle in MeV of various nuclei calculated with a 
modern many-body method with input 3-body yEFT potential. Precision estimates 
are made using an upper and lower bound generated from two different regularization 
cutoffs. Accuracy and precision generally decrease in heavier nuclei. Figure from Ref. 
[ 12 ] 


With confidence in the precision of the many-body methods, the implication is that 
this inaccuracy stems from the input Hamiltonian. Furthermore, this method uses 
an HO-basis SRG-evolved Hamiltonian as input, and there is circumstantial evidence 
that the evolution itself is a major source of the inaccuracy [[26] , 

1.2 Discretization and Numerical Error 

Because we examine the precision of certain aspects of our calculation, a basic 
knowledge of numerically solving integrals and integral equations is critical. As an 
example, we will briefly review numerical quadrature to introduce key terms and 
discuss numerical errors that will be critical in later chapters. Computationally, we 
can approximate an integral as a finite sum of the product of a set of weights (often 
called mesh weights), W{, multiplied by a set of functional values, /) = /(x*), evaluated 
at each element in a set of values within the limits of integration (often called mesh 
points): 

/ b N max 

dx f{x)= Y ( w i fi) + t 

i =0 
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( 1 . 1 ) 





£ is an error which can depend on mesh points, the function being integrated, integra¬ 
tion limits a and b and the integration method. Integration routines are a large topic 
by themselves, and we refer the interested reader to Ref. [S7|. With a finite mesh 
size, one can neither have infinite integration limits nor infinitesimal spacing, thus 
the integration error, £, is created from two sources, which are analogous to infrared 
and ultraviolet errors. First, truncating infinite integrals at a finite value introduces a 
truncation error which will be common in integrals over momentum in later sections. 
Second, there is an error associated with the mesh spacing which would exist even 
over finite integration limits but typically will decrease as our mesh size increases. In¬ 
creasing the range of the mesh but keeping mesh size constant will generally improve 
cut off errors for indefinite integrals, but increase mesh width errors. With these def¬ 
initions, wavefunctions and operators are easily discretized when numerically solving 
integral equations. 

One first must choose a basis to work in; then operators become finite matrices, 
and wavefunctions become vectors; all of which introduce numerical errors related to 
the integration errors previously mentioned. For our calculations, when possible, we 
simply allow the range of the mesh and mesh size to be sufficiently large that errors 
appearing in the 2-body sector from the discretization are negligible when compared 
to the errors we are studying in the 3-body sector. Specifically, we choose mesh 
points and mesh weights from Gauss-Legendre quadrature, which is typically much 
more accurate than Newton-Cotes (constant mesh spacing) or a simple Ricmann- 
summation (for more on these methods, see Ref. [27]). 
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1.3 A Word on Units 


Throughout this paper, units of the potential may look unfamiliar. We choose for 
simplicity to use natural units, h — m— 1. Units depend on basis and dimensionality 
of the system. With this definition, our units all rnnst be powers of fm. To be 
explicit, as a rule, position is always measured in fm, and with natural units this 
means momentum always in fin -1 . The units of the representation of an operator in 
one basis are not necessarily the same as in another, and we will see that the potential 
has units which depend on the number of dimensions and representation. To find the 
units of the potential in a particular basis, we examine the Schrodinger eigenvalue 
equation in that basis and match units for each term. For three dimensions (3-D), 
observing the Schrodinger eigenvalue equation in coordinate representation, 

(Iv 2 + U(r))T(p)=UT(p), (1.2) 

we see that the potential, V(r), and energy, E, both have units of fm -2 , in order to 
match the units of V 2 . The same is true for HO-basis. 

In momentum representation in 3-D, the Schrodinger eigenvalue equation is 

yUfO + J dp V{p,p'yis{p) = E T(p), (1.3) 

and we see that the potential, V(p,p'), must have units of fm, so that combined with 
the fm -3 from dp 1 , the integral has the same units as E T (p) and ^-\l/(p). 

Similarly, in one dimension (1-D) in coordinate representation and HO-basis the 
potential has units of fm -2 . In momentum representation, we again look to the 
Schrodinger equation, 

y*(p) + f dp V(p,p')^(p') = E tt(p), (1.4) 
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and match units to see that the potential now has units of fm 1 . 

1.4 Organization of Thesis 

First, in chapter [2j we review important aspects of similarity renormalization 
group transformations. We will provide details on the different generators used in the 
following chapters and discuss the flow equations that dictate the transformation. We 
also mention a strategy to evolve the 3-body potential while avoiding complications 
from “spectator 5-functions,” which will enter into the discussion of alternate methods 
not used in this study. 

In chapter [3j we re-examine for the SRG the conclusions of Ref. 0 for Vi owk 
potentials, that the potentials must be phase-shift equivalent up to a certain resolution 
scale but also have consistent, explicit handling of the long-range pion physics [6]. We 
use an inverse scattering separable potential (ISSP) to test if universality in potential 
matrix elements emerges at high energies and without explicit pion-exchange terms. 
The ISSP can reproduce all observables in the two-nucleon problem, and we will 
see explicitly that this is not enough for all low-momentum matrix elements to flow 
towards a universal form at finite cutoff. Also, when creating the ISSP we are free 
to choose a binding energy independent of the phase shifts, thus we can see the 
effect of differences in the binding energy on evolved low-momentum matrix elements. 
Furthermore, with the ISSP we are able to generate potentials with local phase-shift- 
equivalence and universal regions, which implies a local decoupling by the SRG. 

To test the idea that the same explicit long-range treatment is required for flow to a 
universal form, we introduce a second simple potential that is phase-shift equivalent 
at low energies and includes explicit one-pion exchange (OPE). We use the model 


11 


proposed by Navarro Perez et al. j281 j2H], which combines the OPE potential with a 
sum of h-shell potentials. This potential replaces the short-range physics with simple 
terms to be fit to phase shifts, while preserving the long-range force. At the close of 
chapter [3j we will set up the 3-body problem with a necessary simplification of the 
phase-shift fitting procedure. 

In chapter [4j we provide details for the treatment of the 3-body problem. We 
discuss the benefits and drawbacks of momentum representation and HO-basis. We 
discuss HO-basis in great detail in the 2-body 3-D sector and for 1-D two- and three- 
body problems. We introduce the imposed truncation errors and implications for 
transformation between bases. When evolving HO-basis potentials with SRG using 
relative kinetic energy as the generator, we identify evolution errors once the decou¬ 
pling scale becomes too low. Because of the singular nature of the kinetic energy, 
accurate transformation between bases is impossible, and different finite bases repre¬ 
sent a slightly different operator due to the truncated basis space. The differences in 
the generator lead to differences in SRG evolution between HO-basis and momentum 
representation. Much of the formalism in setting up the HO-basis will be reviewed, 
however the implications of the basis truncation errors and especially the SRG errors 
will be of critical importance for the following chapter. 

In chapter [5j we create a simple “test laboratory” in 1-D to simplify the 3-body 
problem. We test that our tools reproduce the results of previous papers [3b[|3TJ and 
then introduce a more ideal potential for our needs. We confirm that our model follows 
the same rules for SRG evolution to universal, low-energy potential matrix elements. 
Then we launch into the 3-body sector and discover that the evolution errors of 2- 
body potentials in HO-basis appear at much higher A. These errors prevent us from 
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studying universality of matrix elements in the three-body sector. We discuss the 
origin of the 3-body matrix elements that produce this evolution error. 

In chapter [6j we review our findings and provide ideas for further study. Specifi¬ 
cally, a number of interesting developments are underway in the low-energy nuclear 
theory community that will be able to reconcile the differences in HO-basis and mo¬ 
mentum representation calculations. 

For a list of important abbreviations used throughout, see appendix [Cj 
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Chapter 2: SIMILARITY RENORMALIZATION GROUP 


The similarity renormalization group [111 32] is a continuous series of infinitesimal 
unitary transformations acting on the Hamiltonian. The simplest SRG transforma¬ 
tions can be expressed in differential form as a flow equation: 


dH s 

ds 




( 2 . 1 ) 


where s is a flow parameter mmm- For most nuclear applications to date, the 
operator G s is chosen to be the kinetic energy operator, denoted T. (We will refer 
to G s in this work as the SRG “generator”.) The most commonly used diagonalizing 
generator for non-nuclear applications is known as the Wegner generator p33j. ft uses 
the diagonal of the Hamiltonian defined in momentum representation, H d instead of 
T for G s . Flows using the Wegner generator are indistinguishable from T for the 
range of evolution in the present study but can differ drastically if the SRG cutoff 
becomes very low |22j or if a large-cutoff chiral potential is used [24] , 

The goal of the SRG is to decouple high-energy from low-energy degrees of freedom 
in the Hamiltonian by driving far off-diagonal matrix elements to zero. Instead of s, 
we usually refer to the decoupling scale, A = s _ 2 for T and H d , where A is chosen 
to have the same units as momentum. In the SRG flow with the T generator, the 


dominant term of Eq. (2.1) for far off-diagonal matrix elements is the term linear in 
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Figure 2.1: Matrix elements V(k, k') plotted against kinetic energy, fc 2 , of a realistic 
potential evolved with G s = T re i to different A (oo, 4 fm *, 3 fm 1 , and 2 fm 1 ). The 
matrix elements are scaled by the largest absolute value, then multiplied by 10 to 
emphasize the shape of the SRG-imposed form factor. The colorbar range the first 
scaling, but not the multiplication by 10. The width of the bands is roughly A 2 as 
implied by Eq. (|2.2 ). 


the potential, [[T, V S ],T], where V s = H s — T. If we keep just this term, the flow 
equation is immediately solved for these matrix elements, yielding (with mass m = 1) 


V,(k,k')~V,.„(k,k') e-fr^PP. 


( 2 . 2 ) 


Thus A 2 is roughly the maximum difference between kinetic energies of nonzero matrix 
elements. Once the Hamiltonian is sufficiently evolved to exhibit decoupling, low- 
energy observables can be obtained from a truncated Hamiltonian D33 or one finds 
naturally that a smaller expansion basis is needed for a desired degree of convergence. 


Fig. 2.1 shows potential matrix elements evolved to different A. In all contour plots 
following, we will always normalize the colors to the maximum absolute value of the 
matrix elements (Vjj/max[Vij]), and when noted we further scale the matrix elements 
by a constant for better visualization of shapes. We can clearly see the decoupling 
of energies larger than the width of the red band. When evolving to low enough 
A, the flow equation term is quadratic in V and numerical inaccuracies cause some 
oscillation about the width of the SRG form factor of the dominant term. 
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A nondiagonalizing alternative for G(s) is the momentum-representation block 
generator [23], H bd , defined as: 


H bd = 

PH S P + QH S Q, 

(2.3) 

(p\P\p') = 

"a. 

1 

la 

1 

(2.4) 

Q = 

1 ~P, 

(2.5) 


where 0 denotes the Heaviside step function. Versions of this generator with smoother 
cutoffs exist as well. H b s d matrix elements are the block diagonal elements of the 
evolved Hamiltonian H s , separated at a fixed chosen cutoff parameter A. (That is, the 
generator H bd in a momentum basis is obtained from H s (k, k') by setting to zero the 
matrix elements where k < A and k' > A or k > A and k' < A.) This yields the same 
basic pattern of decoupling achieved with VJowfc Lee-Suziki transformations [611281 15]. 
In fact, the VJowfc an d SRG block diagonal transformations have been shown to result 
in very similar Hamiltonians for the lower energy block if the SRG transformation 
is run to A A [23j. For H bd , A = and represents the maximum difference in 


energy for coupling between the blocks above and below A. Fig. A2 shows a realistic 
potential evolved to different A. Notice that there are two distinct blocks of nonzero 
matrix elements, at k , k' both above A or both below. 

We will see in chapter [3] how different potentials evolve to a low-energy universal 
form. It has been shown that the SRG with the T generator na and Vowfc PEE1E] 
each drive realistic potentials to separate low-energy universal forms, and we will 
show that H bd also drives potential matrix elements to a different universal form in 
the 2-body sector. 
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Figure 2.2: Matrix elements V(k, k') of a realistic potential evolved with G s = H bd 
to different A (oo, 4 fm -1 , 3 fm -1 , and 2 fm -1 ). All A are 1.5 fm -1 , which allows for 
some non-zero matrix elements around the connected corners of the two blocks. 


When evolving in a many-body sector, one can simply plug in the many-body T re j 


and H s into the flow equation, Eq. (2.1). In a momentum representation, 5-functions 
over spectator particles make this strategy numerically impossible, but evolution in a 
discrete basis, such as HO-basis, has no singular 5-functions, thus it apparently solves 
this problem. This is historically the evolution strategy for 3-nucleon potentials and 
will be the strategy we use for our simple system of 1-D spinless bosons in chapter [5] 
Large numerical inaccuracies with this method occur, however, as we will see in 
chapter [4] such that evolution in different finite bases can be different. It is possible 
that the inaccuracies of the Hamiltonian observed in heavy-nuclei problems ra is 
partially due to SRG evolution errors in HO-basis. 

Recent work with SRG evolution in momentum representation takes a different 
approach, explicitly subtracting the singular parts from the flow equation at each 
step [SI CEH1 El]- We can define our 2- and 3-body potential and kinetic energy 
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operators as: 


H — T\ 2 + T 3 + V \2 + W3 + V23 + V123, ( 2 . 6 ) 


Ti2 + T 3 — Ti 3 + T 2 — T 23 + Tf, 


( 2 . 7 ) 


where 12 denotes the relative kinetic energy or 2-body potential between particles 1 
and 2, T 3 is the kinetic energy of the third particle relative to the center of mass of 
the pair, and V 123 is the 3-body potential. The flow equation then becomes, 

dV s dV12 dV l3 dV 2 3 dV 123 _ nr tt 1 o\ 

With some algebra, we can subtract out the 2-body flow equation for each pair, 


dV< 




ds 


— [[Tij , Vij], Tij + Vij]. 


(2.9) 


This subtracts away the “disconnected diagrams” and thus the dangerous spectator 
deltas mm, and we are left with: 

dV 123 
ds 


[P12, V12], V13 + V23 + \ 123] 

(2.10) 

+ [P13, V13]) Vl2 + ^23 + V123] 

(2.11) 

+ [P23, V23], V12 + V13 + V123] 

(2.12) 

+ [[T r euV 12 3 lH s ]. 

( 2 . 13 ) 


This form of the flow equation is not nearly as convenient to implement, as it re¬ 
quires precalculating the 2-body evolved potential with the same grid and weights 
and inserting it at every step of solving the 3-body differential equation. The benefit, 
however is that the 2-body potential in the 2-body space doesn’t have pathologies 
involved with embedding it into a 3-body space (such as spectator 5- functions). A 
very useful feature of this form of the flow equation is that we can see explicitly that 
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the 3-body evolved potential has terms that depend only on 2-body potentials. Thus, 
in a 3-body evolution without an initial 3-body potential, one will be induced. 
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Chapter 3: REALISTIC 2-BODY UNIVERSALITY 


We discuss universality in matrix elements of modern realistic potentials in Sec¬ 


tion |3.1| In Section |3.2[ we provide a working description of the ISSP formalism, 
examine universality in ISSP’s, and discuss the resulting insight into the prerequi¬ 


sites for universality. Section [3731 gives a description of the 5-shell plus OPE potential 
and examines the SRG flow of this potential to a universal form. We also comment 


on the SRG flow of the JISP16 potential. Section 3.5 details the relationship between 
phase-shift and eigenvalue equivalence, which will simplify the 5-shcll fitting proce¬ 


dure for different bases and more bodies. Finally, we conclude in Section 3.4 with a 
summary and the outlook for the three-body problem. 


3.1 Modern Realistic Potentials 


We have chosen a representative phenomenological potential and a set of yEFT 
potentials to evolve and examine in various partial waves. The phenomenological 
potential is Argonne vis (AV18), which employs basis operators in position repre¬ 
sentation and fits the coupling constants to elastic scattering data M. We use 
the N 3 LO yEFT potential from Entem and Machleidt with a cutoff of 500 MeV [35] 
and then five N 3 LO yEFT potentials with various cutoffs from Epelbaum et al. [36]. 
These yEFT potentials have different regularization and phase shift fitting schemes, 
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which creates significant differences in the matrix elements of the potentials. As can 


be seen in Fig. 3.1 all of these potentials reproduce the same low-energy phase shifts. 
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Figure 3.1: Phase shifts of various realistic potentials (see text) in the (a) 1 So, (b) 
3 Si, and (c) 1 Pi partial waves. The shaded regions show the range between the largest 
and smallest phase shifts. The vertical bands indicate the region where phase-shift 


equivalence between the potentials ends, as defined by Eqs. (3.1) and (3.2) 


From Fig. 3.2 one can see that the diagonals of the initial potentials in momentum 
representation are quite different (the differences are particularly evident in lower 
partial waves, so we focus on those). In making these comparisons, we do not single 
out individual potentials but use a shaded region to highlight the range of matrix 
element variation. As advertised, after evolution the matrix elements collapse at low 
momentum to a universal dependence on momentum (the result at fixed A = 1.5 fm -1 


is shown in Fig. 3.3). This feature is not restricted to the diagonal elements; low- 


energy off-diagonal matrix elements of the potentials also evolve to universal values 


(see Fig. 3.15 below). At higher momentum, the potential matrix elements deviate. 

Following Ref. [5] , we compare phase shift and matrix element deviations to iden¬ 
tify the correlations between phase-shift equivalence and matrix element universality. 
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Figure 3.2: Diagonal matrix elements V (k, k ) of various unevolved realistic potentials 
(see text) in the (a) 1 Sq, (b) 3 Si, and (c) 3 Pi partial waves. The shaded regions show 


the range of values and the vertical bands are from Fig. 3.1 
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Figure 3.3: Diagonal matrix elements of various realistic potentials in the (a) 1 So, 
(b) 3 Si, and (c) 3 Pi partial waves evolved by the SRG to A = 1.5 fm -1 . The shaded 


regions show the range of values and the vertical bands are from Fig. 3.1 
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Figure 3.4: Diagonal matrix elements of various realistic potentials in the 1 So partial 
wave evolved by the SRG to A = (a) 5.0 fm -1 , (b) 3.0 fm -1 , (c) 2.0 fm -1 , (d) 1.5 fm -1 
(marked by the vertical dashed line). The shaded regions show the range of values 
and the vertical bands are from Fig. 


3.1 
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Figure 3.5: The spread of diagonal matrix elements of various %EFT potentials (see 
text) in the 1 So partial wave are shown as shaded regions for the unevolved potential 
and then after evolution to A t = 2.5, 2.0, 1.5, 1.0, and 0.5 fm -1 with the T generator 
(red or light gray). These are compared to the spread of the corresponding matrix 
elements for the H bd generator with A = 2.5, 2.0, 1.5, 1.0, and 0.5 fm _1 fm _1 , all 


and the vertical dashed lines mark A t or A. 


evolved to A = 0.5 fm 1 (blue or medium gray). The vertical bands are from Fig. 3.1 
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In Fig. 3.1, we have identified vertical bands within which the phase-shift equivalence 
among the various potentials ends and significant deviation begins. While identifying 
an exact point marking this deviation will be somewhat arbitrary, we can roughly 
choose a normalized width description that is consistent with visual assessments of 
the phase shift plots. In particular, for each partial wave, the vertical band represents 
the region characterized by: 

0.03 < e{k) < 0.1 , (3.1) 


where 


^ _ $high ij' ) t)ow (A) 

e[k> = A 


(3.2) 


The numerator is the range of phase shifts at a fixed k while A is the range of phase 
shifts for the entire universality region. Our studies imply that the precise definition 
of e is not important; as long as it consistently identifies the regions where phase-shift 
equivalence ends it can be used to consistently compare to the regions where the 
universality of matrix elements end. 

Comparing Figs. |3.2 and 3A, we see that while diagonal matrix elements of the 
initial potentials differ significantly in the region where phase-shift equivalence ends, 
this same region corresponds to where the matrix elements have collapsed to universal 
values by A = 1.5 fm” 1 . This suggests the hypothesis that a prerequisite for matrix 
element universality is phase-shift equivalence. Namely, if there are local regions in 
energy in which potentials are not phase-shift equivalent, then there is no universality 


in those regions (this is tested further in Section 3.2). Examining the diagonals of 


the potentials more closely, we observe that for the and 3 Si channels, the lowest 
matrix elements are not exactly the same. This may be a consequence of not evolving 
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A further. From the T generator curves in Fig. 3.5, we can see that the slight width 
of the band decreases as we evolve chiral potentials to A = 0.5 fm -1 . Also, as we will 
see below, differences in the binding energy of the deuteron play an important role in 
the low-energy matrix elements of the 3 Si potential. 


How low must A be before we see universality? Figure 3 A shows the diagonals of 
the 1 Sq potential evolved to four different A values. The vertical bands correspond to 


the same region where phase-shift equivalence ends for the 1 Sq channel as in Fig. 3.1 


while the vertical dashed line shows the value of A. We see in this partial wave 
(and in others not shown as well) that universality in the matrix elements does not 
occur until A approaches the vertical band. A natural hypothesis is that the matrix 
elements will not fully collapse to universal form until A reaches the region of phase- 
shift equivalence. There may be an intrinsic low-energy scale common to each of these 
potentials that determines at what A universality in potential matrix elements will 
appear. A possibility is that this scale is a consequence of explicit treatment of pion 
physics in each of the modern realistic potentials. To test the latter explanation, a 
potential with phase-shift equivalence at much higher momenta and no explicit pion 
physics is required, which we consider in the next section. 

As described earlier, the block-diagonalizing generator will drive the potential 
matrix elements to a different universal form than T. This is illustrated in Fig. |3.5 


with a set of yEFT potentials in the 1 So channel. When evolved to A < 2fm -1 with 
the T generator, the universal form of diagonal potential matrix elements emerges 
over the full region of phase-shift equivalence. For the block diagonal generator with 
A < 2fm _1 , however, only diagonal matrix elements below A become universal and 
with a different flow than the matrix elements evolved with T. The universality is only 
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up to A because this SRG only decouples one block from the other, so matrix elements 
at momenta above A still couple to matrix elements in phase-shift-inequivalent regions 
and therefore do not collapse to a universal form. (Note that in the Vjowfc RG, the 
higher block is set to zero.) We will discuss only G s = T in the rest of this study but 
emphasize that the ideas about universality apply to both generators, although only 
in the low-momentum block for the SRG. 

The region of phase-shift equivalence for the realistic potentials is limited by the 
energies to which they can be fit to clastic scattering phase shifts. The potentials 
are only fit up to the inelastic threshold, about 350 MeV, where contributions from 
pion production become non-negligible. Because of this, if we wish to investigate 
different regions of universality, we must use a method that can ‘fit’ the phase shifts 
in a controlled range of energies. One of the simplest approaches is solving the inverse 
scattering problem with a separable potential, which we consider in the next section. 

3.2 Separable Inverse Scattering Potential 

Instead of fitting coupling constants for predetermined operators to the phase 
shifts, an inverse scattering procedure constructs a potential directly from the phase 
shifts. Separability is just a constraint to define a unique potential, chosen here due to 
its simplicity. For instance, when solving the Lippmann-Schwinger equation, a separa¬ 
ble potential reduces the problem of solving an integral equation to simply evaluating 
an integral. The three-body Faddeev equations also simplify for a separable potential, 
as one of the integrals over internal momenta becomes trivial. A key feature of the 
ISSP for this study is that the potential is entirely created from the phase shifts and 
binding energy of the deuteron; no explicit pion exchange or other physics is imposed. 
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This allows us to determine whether or not universality requires extra physics, such 
as explicit long-range pion terms or other phenomenological considerations. We start 
with a brief summary of the inverse scattering separable potential for two nucleons. 

3.2.1 Formalism 


The form of a rank-n separable potential is: 

n —1 

v = Wi)Aij(vj\ • (3-3) 

i,j =o 

For our purposes a rank-1 separable potential will be sufficient, but future studies may 
benefit from a higher-rank potential. A rank-1 potential in momentum representation 
takes the form: 

V(k,k') = au(k)u(k') , (3.4) 


where a is simply ±1. Details of the rank-1 separable inverse scattering problem are 
well documented mm-, here, we simply state the main results, some limitations, 
and how to work around the limitations (see appendix |A| for derivation). The solution 
to the separable inverse scattering problem is 137]: 


al ,2, k) = ~ k l sin(W) e -A(k) 

k 2 k 


A (k) = -P 


P0 ° dk'5(k' 


7r 


k! — k 


E h = 


h 2 k 2 

2m 


(3.5) 

(3.6) 

(3.7) 


where kb is zero if there is no bound state and equal to the binding momentum for a 
single bound state with binding energy Eb (for a rank-1 separable potential there can 
be at most one bound state, which is the case for the two-nucleon problem in nuclear 
physics). 
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Once v{k) is determined, the entire potential is known from Eq. (3.4). The binding 


energy E & can be tuned independently of the phase shifts. A limitation of rank-1 
separable potentials is that if the phase shift as a function of momentum crosses 
zero, then so too must the potential, and a rank-1 ISSP as defined thus far can 


never change signs if v is real. This point is clear from Eq. (3.8), which follows 
from the Lippmann-Schwinger equation for a separable potential (with standing wave 
boundary conditions): 


-tan (Si(k)) = 


Vi(k, k) 


1,2 dp p 2 Vi(p,p) 
1 + - 7 J p 2 -k 2 


(3.8) 


A zero-crossing in 5(k) corresponds to a zero-crossing on the right side of this equation, 
which can only be achieved by the numerator crossing zero if the denominator remains 
finite. 

Because some of the phase shifts for nucleon-nucleon partial waves exhibit zero 
crossings, we need an inverse scattering potential that allows this feature. We can 
still use the same rank-1 formalism, however, if we split the problem into two energy 
regimes, above and below the zero crossing [38]. Then we can define: 


S<(k) = 

5(k)6(k 0 - k ) , 

(3.9) 

5>(k) = 

8(k)9(k - ko) , 

(3.10) 

V(k, k!) = 

E<(M') + ^>(M'), 

(3.11) 


and determine V < and V > separately using the rank-1 formalism with S < and 5 > as 
input, respectively. We have confirmed numerically that potentials created with this 
prescription accurately reproduce the input phase shifts. 

The method described thus far works directly for uncoupled channels, but for NN 
scattering we must also account for coupled channels, where some further formalism 
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Figure 3.6: Phase shifts using the AV18 potential and the ISSP up to high lab 
momentum k^y, in the (a) 1 Sq, (b) 3 Si, and (c) x Pi partial waves. 
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Figure 3.7: Diagonal matrix elements of the AV18 potential and the ISSP up to high 
lab momentum ky^ in the (a) 1 So, (b) 3 Si, and (c) x Pi partial waves. Cutoff A is in 
units of fm- 1 . 
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is required. For this purpose, we use the Blatt-Beidenharn (BB) convention for phase 
shifts in the coupled channel for our calculations [35J SB]. (In the plots we employ the 
more typically used Stapp-iV convention for the phase shifts for visualization [TU].) 
The BB convention can be summarized as: 

S(fc) = U*{k)A(k)U(k) , (3.12) 

^ / „2i8o(k) q \ 

= (^ 0 e 2i6i(k)J > (3.13) 

Tj/n _ ( cos(e(fc)) sin (e(k))\ . , 

} ~ V-sin(e(fc)) cos {e(k))J ' j 

Here, S(k) is the scattering matrix which relates incoming and outgoing wavefunc- 
tions BB. with k the momentum corresponding to the interaction energy. Then the 
inverse scattering potential can be written as: 

V (k, k') = U f (fc)V(fc, k')XJ(k') , (3.15) 


where 


V(k,k') = 


Vo (k,k') 0 


(3.16) 


0 Vi(k, k') / 

To proceed, one uses the inverse scattering method for uncoupled channels to find 
Vo(k,k') from 6o(k) and Vi(k,k') from 5i(k). The complete potential is then found 
by a rotation by the mixing parameter, e(k). With this complete separable inverse 
scattering formalism, we can now create a phase-equivalent potential at all energies 
in any given partial-wave channel. 


3.2.2 Universality in separable inverse scattering potentials 

We use phase shifts from Argonne Ui 8 to create the phase-equivalent ISSP. In 


Fig. 3.6, we see that the elastic phase shifts are quantitatively reproduced well above 
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the inelastic threshold. We choose Argonne v is specifically because it has phase shifts 
that extend to this high energy, but any realistic potential could be used for starting 
phase shifts. (Note: for simplicity we treat the problem non-relativistically with only 
elastic scattering because we are interested in testing universality and low-energy 
effects, not to have a realistic description of high-energy physics.) The ISSP’s from 
chiral potentials exhibit similar behavior, except that the internal cutoffs drive matrix 
elements and phase shifts to zero at high energies, which is less useful for the present 
investigation. The accuracy of the 1SSP in reproducing phase shifts can be further 
increased simply by using more grid points and increasing the maximum momentum 
if the phase shifts are nonzero above this momentum. 


Figure 3/7 shows the diagonal matrix elements of Argonne V\s and the 1SSP for 
three different partial waves before and after SRG evolution. We observe that after 
SRG evolution to A = 1.5 fm -1 , universality in the diagonal matrix elements also 
extends to the full range of energies. In fact, the only discernible difference in the 
evolved potential diagonals is below the SRG cutoff. Above A the matrix elements in 
the region shown are completely collapsed to universal values. 

Because the binding energy in the 1SSP formalism is independently tuned from 
the phase shifts, we can investigate in the deuteron 3 Si- 3 D! coupled channel how 
universality in potential matrix elements is affected by differences in the bound-state 
energy. Figures [3/8] and [T9] show the effect of phase-shift-equivalent potentials having 
the wrong binding energy. In Fig. |3.8[ the 1SSP is created from the phase shifts of the 
Argonne uis potential in the deuteron channel, but with a binding energy of 0 MeV 
instead of 2.224 MeV. ft is evident that the effect on diagonal matrix elements is 
substantial. The low-energy matrix elements of the bare 1SSP tend towards zero as 
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Figure 3.8: Initial and evolved diagonal matrix elements in the 3 Si channel for AV18 
and an ISSP with a binding energy of 0 MeV. 


the momentum decreases. As the potentials evolve, the diagonal matrix elements are 
driven to universal values except that the ISSP is constrained by its binding energy 
to approach zero as the momentum approaches zero. 


A similar effect can be seen in Fig. T9 where instead of OMeV as input binding 
energy, the ISSP is created with input binding energy of 5 MeV. The ISSP reproduces 
this energy better than 100 eV. This potential is overbound and its lowest momentum 
matrix elements are forced lower than if it had the physical deuteron binding energy. 
Again, the higher momentum matrix elements flow towards a universal form because 
of phase-shift equivalence. These plots show that phase-shift equivalence is not the 
only prerequisite for universality in the diagonal potential matrix elements, but a 
correct binding energy is also necessary. (That is, we need S-matrix equivalence for 
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Figure 3.9: Initial and evolved diagonal matrix elements in the 3 Si channel for AV18 
and an ISSP with a binding energy of 5 MeV. 


negative energies as well.) This may account for the small deviations in the potentials 


at lowest momenta in Fig. 3.2 The 3 D X partial wave plots of the corresponding 
ISSP potentials with different binding energies are indistinguishable. This effect only 
appears in the 3 Si potentials. It is possible that a virtual bound state in the 1 S 0 partial 
wave has a similar effect on the evolved low-momentum potential matrix elements, 
but the ISSP cannot tune virtual bound states and residues in the same way it 
accommodates bound states, thus we do not investigate this point further. 


Next we turn to off-diagonal matrix elements. Figure |3.10| shows the potential 
matrix elements V(ko,k) for k o = 0.1 fm -1 as a function of k for the ISSP and all 
of the realistic potentials evolved to A = 1.5 fm -1 . We can see that although these 
off-diagonal cuts for the modern potentials agree at A = 1.5 fm - , the ISSP matrix 
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Figure 3.10: (a) Off-diagonal SRG evolved potential matrix elements V(k,ko) with 

k 0 = 0.1 fm 1 and (b) unevolved half-on-shell T matrices T(k, k 0 ; k$). In both figures, 
the thick line is from the ISSP while the bands are various realistic potentials. 


elements do not. By using a diagonalizing SRG transformation (that is, G s = T), 
the off-diagonal potential matrix elements are exponentially suppressed. Because of 
this, it appears that the ISSP approaches a universal form, but unlike the realistic 
potentials, there is no finite A at which the ISSP collapses to universal form. Fig¬ 


ure 


3.10 shows low-energy half-on-shell (HOS) T matrices from each of the unevolved 


realistic potentials and the ISSP. We observe that the realistic potentials, which will 
evolve to a universal form, have essentially the same low-momentum, low-energy HOS 
T-matrix elements, while the ISSP does not. This is consistent with carrying over to 
the SRG the suggestion from Ref. [5] that HOS T-matrix equivalence is required for 
off-diagonal universality in Vjow k RG-evolved matrix elements, much like phase-shift 
equivalence is required for universality of diagonals. We only show the 1 S 0 partial 
waves, but the same pattern holds for all partial waves. Clearly, matching observables 
is not enough to produce fully universal potentials after evolution, and in the next 
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section we will examine if matching observables and also including the same explicit 
one-pion exchange potential will be enough for potentials to evolve to a low-energy 
universal form. 
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Figure 3.11: Effects of low-energy phase-shift difference on universality, (a) phase 
shifts, (b) diagonals of potentials 


As a further test, we created ISSP’s using altered phase shifts in localized regions 
of energy to see if the flow to universal diagonal matrix elements is disturbed only 
locally. We use the 1 Pi channel for clarity. Figure 3.1l[ a) shows the x Pi phase shifts 
for Argonne uis and for an ISSP that is phase-shift equivalent except for a Gaussian 
bump that we impose by hand at low energy. In Fig. 3.11[ b) we see that the potentials 
evolve to the same diagonal values everywhere but at low energy. Another potential 
was constructed by creating low- and high-energy regions of phase-shift equivalence, 
and imposing a Gaussian bump (around = 4.0 fm -1 ) to create a difference in the 
intermediate energy phase shifts, see Fig. 3.12[a). In Fig. 3.12[b) the evolution to 
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Figure 3.12: Effects of intermediate-energy phase-shift difference on universality, (a) 
phase shifts, (b) diagonals of potentials 


common diagonal values again works everywhere except near where the phase shifts 
disagree. 

We conclude from these figures (and other tests not shown) that the SRG evolved 
diagonal potential matrix elements are altered only in a region localized near the 
altered phase shifts. This suggests that an SRG softened potential is locally decoupled 
such that the integral in the Lippmann-Schwinger (LS) equation for the on-shell T 
matrix can be truncated as: 


2 r-k +A 

T l (k,k]k?) = V l (k,k ) + -P dpp 2 

n Jk -A 


X 


k 2 — p 2 


(3.17) 


where the lower limit of the integral is taken to be zero if k — A < 0. In Eq. (3.17), A 


represents the local decoupling scale, which we will set to SRG A. (In fact A appears 
to be a conservative upper bound for A to quantitatively reproduce phase shifts.) 
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Figure 3.13: Phase shifts of Argonne Ui 8 potential and truncated phase shifts of 
evolved potentials with A = 4.0 fm -1 . Cutoff A is in units of fm -1 . 
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Figure 3.14: Phase shifts of Argonne Ui 8 potential and truncated phase shifts of 
evolved potentials with A = 1.5 fm -1 . Cutoff A is in units of fm -1 . 
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Figure 3.13 shows phase shifts calculated from Eq. (|3.17|) with A = 4.0 frn 1 in the 


1 S 0 channel for the Argonne Ui 8 potential evolved to three different SRG A’s. These 
are compared to the actual phase shifts of the unevolved potential. We see that with 
this large value of A, the truncated phase shifts for even the unevolved potential are 
largely reproduced and the low-momentum phase shifts from evolved potentials are 
indistinguishable from the actual phase shifts. (The periodicity at high momentum 
for A = 1.5 fm -1 is a numerical grid artifact.) In Fig. 


3.14 


we more severely truncate 


the integral in the LS equation to A = 1.5 fm -1 . We see clearly that the potential 
evolved to A = 4.0 fm -1 is not decoupled enough to reproduce the original phase shifts, 
but the potential evolved to A = 1.5 fm -1 has phase shifts identical to the previous 
plot. This suggests that evolution with T does locally decouple energy scales. 



k [fm 1 ] 



Figure 3.15: (a) Off-diagonal SRG evolved potential matrix elements, (b) Unevolved 
half-on-shell T matrices. In both figures, the thick line is the d-shcll plus OPE po¬ 
tential while the bands are from realistic modern potentials. 
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Figure 3.16: Off-diagonal matrix elements of (a) chiral N 3 LO potentials and the 


JISP16 potential (blue line) in the 1 So channel and the same potentials evolved by 
the SRG to (b) A = 2.0fm -1 and (c) A = 1.5 fm -1 . 

3.3 OPE plus (5-shell 

Here we further test the suggestion that explicit treatment of the longest-ranged 
physics is a requirement for potentials to evolve to a universal form [5]. In particular, 
we develop a simple test potential that is (approximately) phase-shift equivalent in 
the same momentum regions as the realistic potentials but also has the same explicit 
long-range forces. We use the model from Navarro Perez et al. that combines the 
one-pion exchange (OPE) potential with a sum of N h-shcll potentials [28l [29 | in 
each partial wave: 


N 



(3.18) 


2=1 
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The explicit form of the OPE potential is [42], 


V l 7 T (r) 
T(r) 
Y(r) 


3 

rri'i 

7T 

12tt 



2 

ri • r 2 [T(r) 5i2 + Y{r) oq ■ cr 2 ], 


e ~m n r 

m n r 



3 

m n r 


+ 



> 


e 


—m^r 


m^r ’ 


(3.19) 

(3.20) 

(3.21) 


S 12 = 3 (cr x • r)(cr 2 • r) - oq • <r 2 . 


(3.22) 


We choose the {r^} as short-range lengths (under 2 fin), and fit the {g\} to match 
low-momentum phase shifts. For efficiency in momentum representation, we choose a 
different regulator than Ref. | 28l [29] , instead regulating the potential in momentum 
representation with a separable form factor: 

/ reg (M') = e- (fc/A)4 e- (fc ' /A)4 , (3.23) 


for which we choose A = 3 fm -1 . We now have a potential with explicit long-range 
pion terms and adjustable short range terms, which is phase-shift equivalent at low 
momentum to the realistic potentials. 


3.3.1 Universality in OPE plus d-shell 


We can see from Fig.|3.15 that the OPE plus 5-shcll off-diagonal potential elements 


evolve to the same universal form as the modern realistic potentials. Also, Fig. 3.15 
shows the corresponding unevolved HOS T matrices. We see that the OPE plus 5- 
shcll potential has the same low-energy low-momentum HOS T matrix and shows 
a corresponding low-momentum universality in off-diagonal matrix elements. This 
behavior is not unique to the 1 S 0 partial wave, but appears for all partial waves. This 
simple potential explicitly contains only the longest range OPE potential and has 
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very simple short-range terms, but it collapses to the same universal low-momentum 
potential after SRG evolution. Combined with the ISSP results, this is strong evidence 
that the same explicit inclusion of the longest-range contributions to the potential, 
which is reflected in low-energy HOS T-matrix equivalence, is required for collapse to 
a universal form. 


3.3.2 JISP potential 


In principle, a good test of our observations about universality is the JISP16 
potential, which is a realistic potential constructed using the J-matrix version of 
inverse scattering theory [151 i ll] . Because there is no explicit incorporation of a pion- 
exchange tail in the functional form of the potential, we might expect the Hamiltonian 
to exhibit non-universal evolution with the SRG for off-diagonal matrix elements. In 
fact, the unevolved JISP potential is already soft and changes only slightly under SRG 


evolution. But as shown in Fig. 3.16 in the x So channel for a set of off-diagonal matrix 
elements (and true for the diagonal and other partial waves), JISP16 is already close to 
the universal form reached by the chiral N 3 LO potentials. There are still differences, 
but they are small. However the JISP HOS T matrix is also close to the others 
(perhaps as the result of additional adjustments of the potential using the freedom of 
the inverse scattering framework fSj), so there is no inconsistency with our general 
conclusions. 


3.4 Recap and moving into 3-body 


Modern realistic two-nucleon potentials exhibit a flow to universal potential ma¬ 
trix elements under the similarity RG. High and low momenta are decoupled in this 
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universal matrix, allowing us to truncate the matrix and drastically simplify low- 
energy bound state and reaction calculations. Any initial interaction that yields this 
universal matrix after SRG evolution is equally effective. This is of little practical 
importance for the two-nucleon potential, but it could be extremely useful if rnany- 
nucleon potentials display this same type of universality. Producing accurate realistic 
few-nucleon potentials is extremely difficult, and much effort is spent working on pro¬ 
ducing momentum representation matrix elements for additional terms in yEFT. Our 
results suggest that any convenient potential that includes long-range pion exchange 
interactions can be used to produce universal many-nucleon interactions when evolved 
with an SRG transformation. 

Our study of universality for two-body potentials yields the following observations: 

• Inverse scattering separable potentials, with no explicit consideration of long- 
range pion exchange, exhibit a universal collapse of diagonal matrix elements 
after evolution in regions of phase-shift equivalence. 

• If an intermediate region of phase-shift inequivalence is imposed, the collapse 
does not occur in this region, but still occurs in every region of phase-shift 
equivalence. This implies that SRG softened potentials are actually locally 
decoupled in energy/momentum. 

• An incorrect binding energy has a strong effect on the lowest potential matrix 
elements and will prevent flow towards a universal form. 

• While phase-shift equivalence and correct binding energies (i.e., S-matrix equiv¬ 
alence) are apparently requirements for universality in two-body potential ma¬ 
trix elements, the ISSP example shows that these are not sufficient to guarantee 
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a potential that will flow to the same off-diagonal values as conventional realistic 
potentials. 

• However, a potential that reproduces low-energy observables and contains ex¬ 
plicit long-range (OPE) terms does flow to universal form, which is consistent 
with observations made for Viowfc evolution in Ref. [5]. 

• To the extent that low-energy HOS T-matrix equivalence indicates long-range 
equivalence of potentials, it signals off-diagonal universality in evolved potential 
matrix elements. 

• For universality to appear, the SRG decoupling parameter must be sufficiently 
low that potential matrix elements in the low-momentum region of HOS T- 
matrix equivalence are decoupled from high-momentum matrix elements. 

These considerations address the onset of universality for the two-body part of the 
inter-nucleon potential but for a complete discussion we have to consider the full 
many-body Hamiltonian. It is well established that evolution of A induces many- 
body forces of increasing importance [3T1 6) '13] and the SRG transformations will 
only be approximately unitary if they are omitted. This entails a lower limit to the 
region of universality in practical applications. The following chapters will detail the 
framework in which we search for universality in a simple 3-body problem, but before 
that, the next section will detail an important simplification of the fitting procedure. 

3.5 Phase-shift equivalence and eigenvalue equivalence 

In the 3-body problem, the same prescription of fitting h-shcll potentials to phase 
shifts is made much more complicated, because 3-body scattering is much more 
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complicated. It is well known that a relationship exists between the phase shifts 
and the bound states of a system given certain boundary conditions (eg. Luscher’s 
method jlS])- Using this as a guideline, we can confirm that fitting the low-energy 
eigenvalues (including positive eigenvalues) for a Hamiltonian matrix also fits the 
phase shifts in the same energy-regime. Because the eigenvalues are trivial to find, 
many-body scattering formalism is unnecessary for our 3-body study. 

The bound-state eigenvalues, in general, will be basis and mesh independent as 
long as transformations are accurate and numerical precision adequate. This is not so 
for positive eigenvalues, however, which are highly dependent on the basis and mesh 
we choose to work with. In momentum representation, for instance, the scattering 
eigenvalues of the Hamiltonian matrix are almost completely dominated by the eigen¬ 
values of the kinetic energy matrix at large momenta. For this reason, it is convenient 
for numerical precision and visualization to use the difference of the eigenvalues of 
the Hamiltonian matrix minus the eigenvalues of the relative kinetic energy matrix, 
all divided by the mesh weights. We will call this vector, £. 


& = 


eig(H)i - eig(T)i 


Wi 


(3.24) 


£ is clearly defined in uncoupled channels, but how to divide out the weights in 
coupled channels is a complication to this method. Our simple three-body model will 
not have coupled channels, thus we do not develop this method further to take into 
account coupled channels. 

shows £ calculated for the same chiral and phenomenological potentials 


Fig. 


3.17 


as before. Comparing to Fig. 3.1, the same qualitative behavior exists between the 
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Figure 3.17: Hamiltonian matrix eigenvalues minus relative kinetic energy matrix 
eigenvalues all divided by mesh weights (£(&)) of various unevolved realistic potentials 
(see text) in the (a) 1 So, (b) 3 Si, and (c) x Pi partial waves. The shaded regions show 
the range of values and the vertical bands are calculated in the same manner as 
Fig.|3X but with denser meshes for better numerical precision. 


phase shifts and £ (although £ requires more mesh points for precision). We see low- 
energy eigenvalue equivalence in the same regions that exhibit low-energy phase-shift 
equivalence, and equivalence in £ from each potential matrix begins to break down at 
the same energy as the phase-shift equivalence. Therefore, we can safely fit 5-shells 
to low-energy eigenvalues rather than phase shifts. 
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Chapter 4: HARMONIC OSCILLATOR BASIS 


Choosing a convenient basis is a critical step in solving low-energy nuclear prob¬ 
lems. Until now, we have chosen to work with a plane-wave basis with the momenta 
of particles serving as continuous degrees of freedom. Plane waves are of great use in 
simple analytic problems and for uniform systems (e.g. nuclear matter [56]), but for 
finite bound-state calculations, it is far more efficient to expand in a different basis. 
In nuclear few-body bound-state problems, it is most common to use the harmonic 
oscillator (HO) basis. The basis functions of the HO basis are the eigenfunctions of 
the harmonic oscillator Schrodinger equation. 

The HO-basis has a number of useful features which simplify three- and many- 
body problems compared to the plane-wave basis. Most useful to our calculation is 
that the HO-basis is a discrete basis. This simplifies operator equations, which in 
plane-wave basis involve integrals over continuous momenta. Instead of discretizing 
these integrals with momentum meshes and weights, and representing them as matrix 
multiplication, the HO-basis already represents all operators as discrete matrices. 
The HO-basis also simplifies bound states calculations; simply choose your favorite 
eigenvalue routine and input the Hamiltonian matrix; this is true even for many-body 
bound states which require Fadeev formalism in a plane-wave basis An extremely 

difficult aspect of an A-body calculation is the presence of spectator delta-functions 
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that arise from the spectator particle when embedding fewer-body potentials into 


the A-body space (we will discuss this further in section 4.1.3). In the HO-basis, 
the numerically difficult to handle Dirac-delta functions become Kronecker delta- 
functions, but the symmetrization process still connects off-diagonal matrix elements 
in the potential. The simplification of the many-body problem in a HO-basis are 
significant, because the basis size required is greatly reduced. SRG evolution of 
realistic nuclear 3-body potentials in a plane-wave basis is cutting-edge research m 
HU, and these potentials must still be transformed into HO-basis for use in many-body 
methods. 


In this chapter, we start with the formalism for HO-basis expansions in Section 4.1 


We then examine some of the limitations of the finite basis transformation and the 
effects on potential matrix elements in a simple model and with realistic 2-body 
nuclear potentials in Section 4.2| Finally, we will show how the choice of basis affects 


generators and SRG flow of potential matrix elements in Section 4.3 


4.1 Formulation 


Before we jump into our 3-body calculation in 1-D, some of the specifics of HO- 
basis expansions are important to understand. We plan on observing matrix elements 
of evolved potentials, so a firm understanding of exactly how the matrix elements are 
generated is critical. The following section will detail the HO-basis expansion first 
in 1-D for 2 particles, followed by 2 particles in a 3-D partial wave expansion, and 
lastly in 1-D for 3 particles. Because we will examine three identical spin-zero bosons, 
we will build the formalism with that model in mind, but provide some generaliza¬ 
tion to fermions when appropriate. For nonzero spin, both spatially symmetric and 
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antisymmetric states combine with spin and isospin to create total states that are 
antisymmetric (symmetric) for fermions (bosons) under the interchange of particles, 
thus examining spin-zero bosons focuses on a subset of channels used in the more 
complicated nuclear many-body problem. 

4.1.1 HO basis for 1-dimension 2-body 

One can use the HO-basis to expand the wave functions in eigenfunctions of the 
harmonic oscillator Hamiltonian. The Schrodinger equation is 


Hho^u) — E\(j) n ), 

P' 2 1 

Hho = x-h -mCl 2 X 2 . 

2m 2 


It is convenient to define the oscillator parameter (taking m = 1) as 


(4.1) 

(4.2) 




(4.3) 


The solutions of 4.1 in momentum representation are, for integer n an. 


<Pn(P) = 


y/n2 n n\ 


H n (bp)e 


-Ubp ) 2 


(4.4) 


The H n {bk ) appearing in Eq. (4.4) are Hermite polynomials [47]. We see that normal¬ 
izable wavefunctions can be expanded into the orthonormal set of HO-basis functions, 
typically with some truncation, N max , as 


Nrr 


y Cii^i). 


(4.5) 


i=0 


This expansion is exact if N max = oo, but numerically, we must choose some finite 
Nmax which will create truncation errors. We can transform potential matrix elements 
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from plane-waves to HO-basis: 


v nm = (n\V\m), 



(4.6) 


When numerically generating the oscillator wave functions, it is best to use a recur¬ 
rence relationship for Hermite polynomials, 


H n = 2xH n -x + 2(ft - l )// n _ 2 . 


(4.7) 


as multiplying and dividing factorials can quickly lead to numerical round-off errors. 
It is important when expanding with large N max to utilize a finer-grained momentum 
mesh to accurately integrate over the many oscillations of the HO wave functions. 

Once we have the HO-basis potential matrix elements, we can also simply con¬ 
struct the relative kinetic energy matrix using ladder operators. We use the iden¬ 
tity [H], 

( 4 - 8 ) 

in the HO-basis and find the relative kinetic energy (see appendix |B|), 

(n\T (2) \m) = -^((2n + 1 )<5 n , m + \/{m + l)(m + 2) h n , m+2 

+ \/(n + l)(n + 2) 5 n +2,m)- (4-9) 


In (4.9), the superscript (2) means that this is the 2-body relative kinetic energy, 
and Sij is a Kronecker delta. We see that the relative kinetic energy takes an infinite 
tri-diagonal form when the basis is not truncated. As we will explore later, truncating 
T in either momentum representation or in an HO-basis will have consequences when 
transforming between the two bases. 
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For indistinguishable particles, we must symmetrize (antisymmetrize) the Hamil¬ 
tonian for bosons (fermions). This is done quite simply for two particles in the 
HO-basis. The 2-body problem can be split into purely even and odd states which 
are orthogonal to each other. For spin-zero bosons, we keep the symmetric states 
(even-n) and omit all odd states (odd-n), and for spin-), fermions we can mix the spa¬ 
tially symmetric even-n (antisymmetric odd-n) states with the antisymmetric singlet 
(symmetric triplet) spin states to create overall antisymmetric states. 

The usual method for finding bound state energies will be to start with an analytic 
form of the potential in momentum representation, and then transform into HO-basis 


with Eq. (4.6). Then we add the HO-basis relative kinetic energy matrix to form 


the full Hamiltonian matrix. The Hamiltonian is diagonalized with any standard 
routine on a computer. The 2-body problem can be solved numerically in momentum 
representation just as easily as in a HO-basis, but as we will see, adding a third particle 
complicates the problem in momentum representation, and thus the HO-basis greatly 
simplifies the calculation. 

4.1.2 HO basis for 3-dimension 2-body 


Before adding a third particle, it is important to note that the HO-basis is used in 
realistic, 3-D calculations as well. The 1-dimensional formalism can be extended to 3 
dimensions. We expand into partial waves which introduces a new parameter, /, the 
orbital angular momentum. The solutions of the 3-D partial-wave harmonic oscillator 
Schrodinger equation in momentum representation, with energy, E = 2n + l + |, 
are [SSj: 


”0nZm(k) ( l) 


2 n\b 


{kb ) +1 (kb) 2 


r(n + Z+f) k 


-L\f> ({kb) 2 )Y lm {k), (4.10) 
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l_ I_1 

where L n 2 ( x ) are the generalized Laguerre polynomials, T(a;) is the gamma function, 
and Y lm ( k) are the spherical harmonics. 

4.1.3 HO basis for 1-dimension 3-body 

The three-body problem introduces a great deal of complexity into solving bouncl- 
and scattering-state problems. Here we examine only bound states of our 1-D sys¬ 
tems. The first step in simplifying the problem, much like in two-body systems, is 
changing the degrees of freedom from the individual particles’ positions or momenta 
to a reduced set of degrees of freedom by subtracting out the center-of-mass motion. 
Taking the position of each particle as xi,x 2 ,x 3 and assuming equal-mass particles, 
we can define a set of Jacobi coordinates (r, R)i, where i indicates the spectator par¬ 
ticle. Unless otherwise noted, we will choose particle 3 as the spectator particle, and 
then we define: 

r = =(x 1 ~x 2 ), (4.11) 

R = \j\ Q(®i + x 2 ) -x^j . (4.12) 

We see that, up to a constant, r corresponds to the separation of the pair of particles 
1 and 2, and that R corresponds to the distance of particle three from the center of 
mass of particles 1 and 2 (if all masses are equal). The choice of the constant is to 
ensure that the relative kinetic energy operator is symmetric in terms of the conjugate 
momenta p and q, defined in terms of single particle momenta, ki, as: 

p = - k 2 ), (4.13) 

q = + ^ ” ks ) ' ( 4 ‘ 14 ) 
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Thus the total relative kinetic energy operator is: 


r > 2 n 2 
J<3) _ R _|_ 

2m 2m ’ 


(4.15) 


where m is the mass of each particle. Figure |4.1 shows the Jacobi coordinates for 
a 2-D arrangement of three particles. We show two dimensions for better visualiza¬ 
tion of each vector. From here we can transform from plane-wave basis in Jacobi 
momenta (basis states | pq) = \p)\q)) to HO-basis (basis states | n p n q ) = \n p )\n q )) by 
transforming each momentum separately, namely: 


( pq\n P n q ) = <f> n P {p)<t> nq (q) (4-16) 

(n p n q \V\n' p n' q ) = j dp dq dp' dq (n p n q \pq)(pq\V\p q)(p'q\n' p n' q ). (4.17) 

It is important to make our N max truncation on total oscillator number, rather than 
on each n v and n q independently. This is equivalent to making a truncation in total 
energy. The truncation is important later to ensure that transformations between 
bases with different spectator particles are unitary, which is required for the sym- 
metrization process. 

Now we can emphasize the difficulty with embedding 2-body potentials into a 3- 
body plane-wave basis, and the simplification of HO-basis. The 3-body Hamiltonian 
has three different 2-body interactions; one for each unique pair of particles. If we 
choose to work in Jacobi momenta with particle 3 as the spectator, the interaction 
between particles 1 and 2 has a simple form: 

V^ 12) {p,q]p',q') = V {2] (p,p')5(q-q'). (4.18) 

Here, the superscript, (2,12) refers to the 2-body pair interaction between particles 
1 and 2 embedded into the 3-body basis. 
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Figure 4.1: Jacobi coordinates, u\ = \/2r and u\ 
vectors in 2-D. 


[3 

V 2 


R, and particle positions as 
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Thus we see the spectator delta-function for the non-interacting spectator particle. 
In this basis, however, the other two pair interactions are more complicated. There 
is a linear transformation between sets of Jacobi momenta with different spectator 
particles, so we can write each Jacobi momentum with a given spectator index as a 
function of both Jacobi momenta with a different spectator index; 

PiiPjQj) = a Pj + POj, (4-19) 

QiiPjQj) = lPj + <%• (4.20) 

Here, the subscript on momenta refers to the spectator particle, and thus which set 
of Jacobi momenta p or q belongs to. 

The values of the constants {a, /3, 7 , J} are calculated given the definition of Jacobi 
momenta. Now, the other two pair-wise interactions are written in their own Jacobi 
momenta, which are taken as functions of the first set of Jacobi momenta: 

V {3 ’ 23 \p,q]p',q') = H (2) (pi(p,q),p[{p,q)) 3 (qi(p,q) - q[(p,q)) , (4.21) 

V {3 ’ 31) (p,q;p',q') = V {2) (p2(p,q),p' 2 (p,q)) 3(q 2 (p,q) - q' 2 (p,q )) • (4.22) 

This often requires interpolation when solving numerically, as the Jacobi momenta for 
different spectators cannot be all set on the same momentum mesh and thus must be 
integrated out [Till 30J. There are other subtleties in constructing the potential matrix 
involving consistent domains for momentum integrals [TU], and after the potential is 
constructed, one must use Fadeev methods; finding poles in the 3-body multi-channel 
T-matrix; to find the bound state [16] . 

In the HO-basis, the delta-functions are not Dirac delta-functions but rather Kro- 
necker delta functions. Thus, the 2-body potential between particles 1 and 2 can be 
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embedded into a 3-body basis with spectator particle 3 as: 


y (2 ' 12) , , 

n p n q n' p n' q 


v 


( 2 ) 


(4.23) 


Because we are interested in spinless bosons for our 3-body study, we must symmetrize 
our basis. We follow the same strategy as Jurgenson in Ref. EU- First, we symmetrize 
the 2-body states. In the 2-body system, the symmetrizer is, 

S^ = ^(l+P 12 ), (4.24) 

where P \ 2 is the permutation operator between particles 1 and 2, and 


P 12 \n p ) = (-1H n p ). (4.25) 

Now, the symmetrized states are the eigenstates of the symmetrizer with eigenvalue 
one. We can clearly see that n-even states fulfill this requirement, thus to symmetrize 
the 2-body sector, we omit n-odd states. We restrict our basis states to \N v n q ) = 
in)in, w h ere N p is an even integer. We can then apply the 3-body symmetrizer, 

S (3) = -(1 + P \ 2 + P 2 3 + P12P23 + P23P12 + P12P23P12), (4.26) 

6 

upon our partially-symmetrized basis states. Because \N P ) is an eigenstate of P \ 2 
with eigenvalue one, the symmetrizer reduces to 

S (3) = ^(1 + 2P 23 ). (4.27) 

Refer to appendix [B| to calculate {N p n q \P 2 3 \N! p n' ( ^, after which the HO-basis matrix 
elements of the 3-body symmetrizer, follow directly. 

It is convenient to rearrange our basis once more, because our truncation is in 
total N 3 — N p + 7 i q . We order states by energy and a degeneracy index, rather 
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than the two oscillator parameters. Doing this reorganization reduces the number 


of matrix elements by omitting many that have value zero in our total-n truncation 
scheme, and it also organizes the matrix elements into energy blocks. This is not a 
transformation, but just a reorganization of matrix elements and a shrinking of the 
matrix by removing zero-valued elements. The new basis can be written as | A^ 3 , z), 
where i is simply a label for the degenerate states at each energy. For our matrices, 
we use i = n q . Finding the set of eigenvectors with eigenvalue one of the 3-body 
symmetrizer in HO-basis will give us the coefficients of fractional parentage, 


(N 3 ,i\N^i') s . 


(4.28) 


The whole process can be set up numerically as a series of matrix multiplications 
to first transform from Jacobi plane-wave basis, \pq), then to the energy-truncated, 
partially-symmetric HO-basis, |jV 3 , i ), and finally to the fully-symmetric energy-truncated 
basis, \N 3 ,i)s- The embedded 2-body potential in the fully-symmetric basis is iden¬ 
tical under any interchange of particles, thus the problem of separately treating each 
pair interaction is solved. One simply must multiply by a combinatoric factor, Q), 
for the total embedded 2-body potential. Any explicit 3-body potential will receive 
the same treatment, transforming from momentum representation to a partially sym¬ 
metric HO-basis, to the full symmetric basis, but explicit 3-body potentials have no 
combinatoric factor. 

We must also build the relative kinetic energy operator in HO-basis with ladder 
operators in both the n p space and the n q space (see appendix [B] for full calculation). 
This matrix reduces to: 



(4.29) 
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We create this matrix in the partially-symmetrized \N 3 ,i) basis, then symmetrize by 
transforming into the symmetric basis, just like for the potential. 

Once the embedded 2-body potentials, the explicit 3-body potential, and the 
relative kinetic energy are all in the symmetric HO-basis, eigenvalues and eigenvectors 
can be found by a simple matrix diagonalization routine. To extend the process to 
fermions, one needs only to use the odd n p states, and then find the eigenvectors of 
the antisymmetrizer with eigenvalue one. The antisymmetrizer is constructed simply 
by multiplying each permutation operator in the symmetrizer by (—1). For fermions, 
one can attach a spin-like variable to emulate spin-| particles, which complicates 
the antisymmetrization process. For instance, in the 2-body sector the spin singlet 
(S = 0) is odd and must be paired with even n p , and the triplet with odd n p . Thus, the 
potential would appear band-diagonal in spin-space. The 3-body antisymmetrization 
operator becomes more complicated, but the process to find antisymmetric states 
remains the same. 

4.2 Basis transformations, truncations, and cutoffs 

We can now examine some subtleties in using HO-basis expansions for few-body 
calculations. As mentioned before, choosing a finite truncation, N max , is required to 
numerically solve problems, but also creates truncation errors. Some of these errors 
can be negated by tuning the HO parameter and keeping basis size sufficiently large, 
but certain operators are only accurately transformed numerically with an infinite 
basis (e.g. singular operators) . 
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4.2.1 Cutoffs imposed by HO-basis expansion 

Because of the popularity of the HO-basis in realistic nuclear calculations, much 
work has been done to understand the basis truncation errors. As proposed in 
Refs. |3H, SO], using a finite HO-basis imposes both an infrared (IR) and ultravio¬ 
let (UV) cutoff. We can think of these cutoffs as a length cutoff and momentum 
cutoff. For 1-D these are, 


Lho by/^N max T 5 , 

A ho — - yj 2 N max + 5 . 

Or, for 3-D S-waves, 

Lho = by/ 4 N max + 7 , 
A ho — \f ^N rnax + 7 . 


(4.30) 

(4.31) 


(4.32) 

(4.33) 


As long as all non-negligible potential matrix elements are at separations below Lho 
in coordinate representation and momenta below A ho hi momentum representation, 
then we can transform between plane-wave and HO-basis with high precision. This 
gives us a method by which we can calculate the requirements on b and N max for 
accuracy in our basis transformations. If our basis size is too small, or b is not 
tuned properly, then the HO-basis truncation will truncate non-negligible matrix 
elements, and thus transformation back into plane-waves will produce oscillatory 


errors in the matrix elements. These effects are shown in Figs. (4.2 4.5). We can 
see an indication of an inaccurate basis transformation in the matrix elements of 
the potential in HO-basis if the matrix elements along the edges are nonzero. This 
implies that we have truncated nonzero matrix elements and thus transformation to 
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momentum representation will produce errors in the potential matrix elements. The 
transformation back into plane waves is missing the high-rt basis functions, and thus 
we see oscillations in the matrix elements. 

To illustrate this point, we choose a very simple potential matrix, 


V(p,p') = e P 2 e ? 2 ; (4.34) 

V(r, r') = e -Z 2 ~e -r 2 ~. (4.35) 


We choose this form because of the ease of Fourier transforming, the presence of a 
cutoff momentum, and the symmetry in momentum and coordinate representation, 
which allows us to illustrate the requirements both on A HO and L ho . Observing 


Fig. |4.2[ we see the off-diagonal matrix elements of the potential in momentum repre¬ 
sentation and a contour plot of the potential in ffO-basis. We choose b = 1.0 because 
of the symmetry between momentum and coordinate representation, and N max = 6 
for convenience (for this simple case, it is clear that only one oscillator is necessary). 
This potential is essentially an outer product of the n — 0 HO wave function, thus 
only n — 0 is required, but taking b close to, but not equal to 1.0 continuously adds 
small values to higher-n matrix elements. 

This choice of b and N max allows for accurate transformation back and forth 
between HO-basis and plane-waves. This is no surprise, as calculation of the HO- 
induced cutoffs yields k HO = L HO = vTf. Observing the off-diagonal slice of the 
potential, this is in the area where potential matrix elements are negligible. Because 
of the symmetry between momentum and coordinate representation for this potential, 
and 6 = 1, the same is true for coordinate representation. 

Now, if instead we transform to an HO-basis with b = 2.0, and the same lV max , we 
can calculate the HO-induced cutoffs and find A ho = and Lho — 2\/l7. The IR 
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V(0,k) 



Figure 4.2: (Left) Off-diagonal matrix elements of our test potential matrix elements 
in momentum representation (solid blue line) and plotted again after transformation 
to HO-basis (b = 1, iV max = 6) and back (dashed red line). (Right) HO-basis potential 
matrix elements (b — 1, iV max =6). To best show occupancy, we plot 14, n' /14iax, the 
matrix elements divided by the max value. 
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error actually improves, but the induced momentum cutoff truncates non-negligible 
matrix elements, thus there will be some error in transforming from HO-basis to 
plane-wave basis and visa versa (actually, there is no error transforming from coordi¬ 
nate representation to HO-basis, but going back, or any transformation to and from 
momentum representation will not be accurate). We can see a contour plot of the new 
HO-basis potential matrix elements and a cut of the momentum representation ma¬ 
trix elements of the original potential and twice transformed potential (transformed 


from momentum representation to HO, and then back) in Fig. 4.3 



We can clearly see that the HO-basis matrix elements are non-negligible at the 
edge of the basis truncation (and will see later that indeed we have truncated nonzero 
elements). The effect of this truncation is clear when observing the cut of momentum 
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representation potential matrix elements; the cutoff forces higher momentum matrix 
elements to zero and oscillatory errors appear at momenta below the cutoff. 


The message is similar if we choose b = 0.5 as well. Now, the UV cutoff is 
higher, thus transformation from momentum representation to HO-basis will actually 
be more accurate, but transformation back, or to/from position representation will 


have truncation errors. We can see this in Fig. 4.4 



The truncation introduces IR errors, and sure enough, observation of the cut of 
momentum representation matrix elements reveals the characteristic oscillation of a 
long-range truncation. A plot in coordinate representation would look very much the 


same as Fig. 4.3 Once again, the HO-basis matrix elements are nonzero at the edges, 
and thus there will be an incomplete transformation. 
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Often, the choice of b is made to minimize UV or IR errors, or the error on an 
observable, so to increase accuracy of a calculation we must increase N max . For our 
test potential, we saw that b = 1.0 is ideal for basis transformations, but if we want 
to get the same UV cutoff and accuracy for b = 2.0, we must quadruple N max . 



Figure 4.5: Same as Fig. 4.2, but for b = 2, N max = 24 


Figure [475] shows the HO-basis potential matrix elements and a cut of the momen¬ 
tum representation matrix elements of the original potential and the doubly trans¬ 
formed potential. We can clearly see that our new N max truncation is large enough 
such that only zero-valued matrix elements are omitted. Increasing accuracy in HO- 
basis transformations by increasing N max is a costly game, as the scaling of the cutoff 
is proportional to i\U ax . 

In practice, when observing potential matrix elements, we can look for the same 
indications of inaccurate basis transformations due to truncation. If we observe 
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unexplained oscillations in momentum representation potential matrix elements, or 
nonzero valued matrix elements at the edges of our HO-basis potential matrix, then 
we know that our truncation has induced errors in the matrix elements. 


Our goal is to observe universality in the matrix elements, thus oscillatory errors 
will make it impossible to follow the same strategy as before, so we require very 
accurate basis transformations. Often, such accuracy in basis transformations is 
unnecessary for calculating certain observables (e.g. ground state energy of a 2- 
body system), but accuracy in all non-zero matrix elements and likely all many-body 
observables requires larger HO-bases. 

For our calculation, our choice of iV max will be large enough such that there are no 
transformation errors identified as before, and b will be such that the 2-body ground 
state binding energy is accurate to several digits (which will mean low 1R errors). 
This will allow us to choose a preliminary basis size to accurately transform between 


plane-wave and HO-basis for unevolved potentials, but as we will see in Section 4.3 


accurate SRG evolution of the potentials requires more careful consideration of oscil¬ 
lator parameter and basis size. 


4.2.2 Transforming singular functions 

Numerical transformation of singular functions, specifically the Dirac 5-function 
between plane-wave and HO-basis is one major limitation of choosing a finite N max 
truncation. Our 5-shcll fitting procedure is introduced using Dirac 5-functions in 
coordinate representation, and the relative kinetic energy operator introduces a 5- 
function in momentum representation. The fitting procedure requires only narrowly 
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peaked functions which will be fit to data, thus we can instead use regulated 5- 
functions and completely bypass basis-transformation complications, but altering the 
5-function in the kinetic energy imposes unavoidable errors. 

To better understand the limitations of singular functions, we examine the 5- 
functions in coordinate representation. One may assume that a simple solution would 
be transforming to HO-basis analytically and then truncating the space. We see this 
is simply done: 

A nm (a) = J dr dr' </> n (r)5(r - a)5(r - r')0 m (r'), (4.36) 

A nm (a) 0 n (ci)0 m (ci). (4.37) 


The trouble lies in the normalization factor in the oscillator wave-functions. If we 
wish to keep errors negligible, then we should choose a sufficiently large N max such 
that the truncated matrix elements are negligible. Using the closed form for Hermite 
polynomials with argument zero (called Hermite numbers) [51J [52], 


Anm(O) 


A nm (0) 


0n(O)0 m (O), 

(n — 1)!! (m — 1)!! 
■\fn\ yfm\ 


(4.38) 

(4.39) 


where nU denotes the double factorial. Utilizing an identity of the double factorial m 
E2] USES, this can be expressed for non-zero, even n and m as 

*-<“>^ (II f?) (II f ?)• 

This is extremely slow to converge and thus any finite basis HO-basis we hope to use 
will have truncated non-negligible matrix elements. 

Instead, we choose to use regulated 5-functions for our fitting procedure. We 
do not require explicit 5-functions for our fitting procedure, simply narrow-peaked 
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functions in coordinate representation, and by regulating, we ensure that we avoid 
the complications of non-zero matrix elements at the edge of our basis truncation. 
We regulate the coordinate representation 5-functions by first Fourier-transforming 
to momentum representation and then imposing a soft momentum cut-off with a 


regulating function; just like before in Section 343 Then, we transform this regulated 
5-function to HO-basis. If the regulating function is not too sharp, and the cut-off 
is sufficiently small, then this regulated 5-function can be transformed between HO- 
basis and plane-wave basis with high accuracy. Because the fitting routine we use 
does not rely on the functions specifically being 5-functions, this regulated 5-function 
works just the same in creating the 2-body potential. 

Alternatively, the relative kinetic energy operator in a plane-wave basis is a 5- 
function in momentum representation, 


(p\T\p') = J^5(p - p'). (4.41) 

We can not accurately replace this operator with a regulated 5-function, so we cannot 
accurately transform the relative kinetic energy matrices between plane-wave and 
HO-basis. Regulating (or simply imposing a momentum cutoff) causes errors in all 
of the liigh-n, HO-basis kinetic energy matrix elements which will make calculations 
for observables impossible. Instead, we use the analytic expression in each basis. 
Because of basis truncation, the operator represented in one basis will not be identical 
to the representation in another basis. Because the kinetic energy is tri-diagonal 
in HO-basis, this method isolates the truncation error to the two matrix elements 
above our truncation which mix the omitted and retained matrix elements. Using 
this method, we are able to accurately reproduce observables in the 2- and 3-body 
sectors. Because of this restriction, we can only transform potentials between bases 
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and not the full Hamiltonian. To transform a Hamiltonian in one basis to another, 
we subtract the relative kinetic energy matrix calculated in the original basis, then 
transform the potential and add back the relative kinetic energy matrix calculated in 
the new basis. In the next section, we will see that differences in T re i in each basis 
will create differences in the SRG evolution of potentials. 

4.3 SRG differences in plane-wave and HO basis 

Because of the subtleties in transforming certain operators from momentum rep¬ 
resentation to HO-basis, certain generators of SRG flow become impractical when 
evolving in HO-basis. For instance the Wegner generator, G(s) = FO(s), is basis- 
dependent, as the diagonal potential matrix elements in momentum representation 
transformed to HO-basis is no longer a diagonal matrix. If we were to use the Wegner 
operator, at each step in the SRG evolution routine, we would have to transform the 
potential to the plane-wave basis, then transform only the diagonal elements back into 
HO-basis to form the generator, making the evolution computationally cumbersome. 
The best candidates for evolution in HO-basis are well defined entirely in HO-basis 
at every A. For this study, we will use T rel , as it is completely defined in HO-basis 
and also is not A dependent, although we have discussed that the basis truncation 
makes T re \ different in different bases. 

We have seen earlier that evolving potentials with T re \ as a generator in momen¬ 
tum representation softens the potential. This has the effect of reducing the required 
N max for accurate transformation to HO-basis. One can truncate the decoupled high- 
momentum matrix elements of a SRG-softened potential, effectively reducing momen¬ 
tum cutoff, Ap 0 t, above which matrix elements are zero. This means that A ho can 
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be smaller, and N max and b can be tuned such that a smaller basis has no additional 
truncation error. 

Because of the basis truncations and finite mesh sizes, the T re \ matrix in plane-wave 
and HO-basis is no longer identically the same operator, and we should expect to see 
differences in evolved potential matrix elements at some point in the SRG evolution. 
In momentum representation, A pot is effectively reduced by SRG evolution, but what 
about in HO-basis? A very simple argument can be used to at least suggest that A pot 
in momentum representation cannot increase with SRG flow with generator T re i (or 
any generator diagonal in the basis). We can define a corresponding N pot < N max 
in HO-basis corresponding to some value of n above which all matrix elements are 
identically zero, but at which the potential is nonzero. The same argument shows 
that under SRG flow with generator T re \ in HO-basis, there are terms which act to 
increase N pot . 

Take four matrices, D tJ , , M t j and V tJ corresponding to a diagonal matrix, tri- 

diagonal matrix, and two different truncated matrices respectively. The truncated 
matrices are constructed so that we can define some n such that 


Vij, Mij = 0 if i or j > n. (4.42) 

Then, looking at what happens when V t] is multiplied by Dij,Tij , or M tJ 


(DV)„ 

di V%j 5 

(4.43) 

(TV),, 

diVij + bjVij —i CjVi—\j, if i, j 1 , 

(4.44) 

(MV),, 

= J2 M “ v v- 

i 

(4.45) 
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The values of coefficients a*, b, n Ci and di are generally not important. We see then 


(DV)( n +l)j 

o" 

III 

i—l 

+ 

(4.46) 

{TV)(n+l)j 

“f - biV(n-\-l)(j—l) CiVnj ^^nji if ^ j 1? 

(4.47) 

(MV)( n +l)j 

= M( n+1 yVij = 0 . 

(4.48) 


Therefore, the only one of these matrix multiplications that increases the value of n 
is TV. The same is true for reversed-order multiplications. If we examine the RHS 


of the flow Eq. (2.1), we see that it is comprised entirely of these three operations 
and their reverse orders. Particularly, the RHS of the flow equation in momentum 
representation has no tri-diagonal terms, thus elements with index above n will never 
change from zero under SRG flow. If we extend this from a discrete matrix back to 
a continuous potential, then, A pot will not grow as we evolve in momentum repre¬ 
sentation with T re i. Alternatively, because T ie \ in HO-basis is tri-diagonal, SRG flow 
can increase N pot at each step. Indeed, when we examine the evolution of potentials 
in HO-basis, we see that they spread to higher n-values than required for unevolved 
potentials. We shall see that SRG evolution in the HO-basis can push the truncation 


requirements higher than we observed for unevolved potentials in Section 42 It is 
also important to remember that the SRG is a continuous set of unitary transforma¬ 
tions, so even if the HO-basis truncation is too small and will create errors in matrix 
elements, it will still reproduce the same eigenvalues as the unevolved potential up 
to numerical errors in the algorithm. 

Another source of differences between evolution in each basis is the finite momen¬ 
tum mesh size. For larger N max , we must choose a momentum grid with enough points 
to accurately integrate over the many oscillations of the oscillator wave function. As 
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the potential is evolved in momentum representation, however, the high-momentum 
matrix elements become sharply peaked on the diagonal. For this reason, more points 
must be included in our momentum mesh for an accurate transformation to HO-basis. 
This manifests as IR-errors, as the momentum difference between mesh points in mo¬ 
mentum representation is related to a long-range cutoff in coordinate representation. 




k 


Figure 4.6: (Left) HO-basis (b = 1.0, lV max = 6) matrix elements of our test potential 
evolved to A = 2.0 in HO-basis and scaled to 10 times their absolute maximum value. 
(Right) Vertical slice of potential matrix elements evolved to A = 2.0 in momentum 
representation. The solid black line is calculated completely in momentum represen¬ 
tation, the dashed red line is evolved in HO-basis, then transformed to momentum 
representation, and the dashed blue line is evolved in momentum representation then 
evolved to HO-basis and back. 


Figure 4.6 shows a contour plot of HO-basis potential matrix elements and a ver¬ 
tical slice of momentum representation matrix elements for our simple test potential 
evolved to A = 2.0. We can see that the HO-basis matrix elements along the edge 
of our truncation are small, but nonzero. Looking at the momentum representation 
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matrix elements, we see very small differences emerging at some momenta. For this 
test potential, SRG evolution in each basis is the same until A = 2.0. If we evolve 


further to A = 0.9, shown in Fig. 4.7 , we see that evolution in each basis becomes dif¬ 
ferent. At A = 0.9, the HO-basis matrix elements have increased along the edge of the 
basis space. The errors in momentum representation potential matrix elements have 
also increased. The dashed blue line, which is produced by evolving in momentum 
representation then transforming to HO-basis and back, shows that there are errors 
from the basis transformation. The dashed red line, which is produced by evolving in 
HO-basis and then transforming to momentum representation, differs from both the 
black and dashed blue line, showing that the evolution is different in each basis. 

Unlike the unevolved imposed cutoffs HO-basis, the new requirements on basis 
truncation and oscillator parameter for accurate SRG evolution are not analytically 
understood, but a larger basis size allows for accurate SRG evolution to lower A. For 
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this reason, some choose to transform between HO-basis with different b and N max 
dynamically during their evolution [53] • For our study, we will simply start with a 
sufficiently large basis and tuned oscillator parameter such that we do not observe the 
oscillatory errors and nonzero edges in the 2-body evolved potential matrix elements. 

4.4 HO-basis Recap 

In this chapter, we have explored the many sources of errors when transforming 
between finite bases. Transforming to HO-basis adds a cutoff in both momentum 
representation (UV cutoff, A ho) and coordinate representation (IR cutoff, Lho )• 
In the 2-body sector, if these cutoffs are large enough such that only zero-value 
matrix elements are omitted, then there will be no truncation errors in transformation 
between HO-basis and plane-wave basis. 

Because singular functions cannot be accurately transformed between bases, we 
choose to use the analytic expression for the kinetic energy in each basis. The kinetic 
energy matrix in a finite HO-basis is not the same operator as in a finite momentum- 
representation plane-wave basis, because of the truncation. For this reason, SRG 
evolution performed in each basis will have differences that increase as A decreases, 
even if the unevolved potentials can be transformed between the basis without trun¬ 
cation error. If we increase basis size, we can reduce the SRG differences in each 
basis, thus we can use the 2-body sector to define an N max large enough such that 
SRG errors are negligible for the A we evolve to. 
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Chapter 5: TOWARD 3-BODY UNIVERSALITY IN A 


“1-D LABORATORY” 


With the formalism in place, we can begin the 3-body calculation in 1-D. We 
will introduce a potential utilized in previous studies j30j [31] to confirm binding 
energies, then introduce a better potential for use in our study. We will show that 
SRG evolution errors appear in the 3-body sector at much higher A than expected 
from the 2-body sector analysis of the previous chapter. We can identify these errors 
in the potential matrix elements as truncation errors originating from the spectator 
(5-functions. Although the potentials reproduce 3-body observables, each step of the 
SRG evolution compounds the truncation errors, causing increasing errors in the 
evolved potential matrix elements as A decreases. These errors make our study of 
universality in 3-body matrix elements impossible, but understanding them brings us 
closer to precision Hamiltonians for many-body nuclear problems. 

5.1 Model Potentials 

First, we will consider as a starting point the potential proposed in Ref. [SDj, 31]. 



(5.1) 


where r is the separation between particles. This potential is easily transformed into 
momentum representation. We must be consistent in our 2- and 3- body bases, so we 
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will use the Jacobi-coordinate, l = r/\J 2 and conjugate momentum, p, instead of the 
usual 2-body center of mass momentum. Then, we see 


V 2 (p,p') = 


Vi(p,p') = 


^ e~ ll{p - p,) V 2 (y/(2)l) 

Z7T 


<*i ( p-p f )) 2 


(5.2) 

(5.3) 


ti^V2 ' 

We can confirm the binding energies of Ref. |30j for different choices of a* and q, 
using b = 0.5 fm and N max = 120. These results also employ a 3-body potential of 
the form, 

2,2 /2 , /2 . 

T r / / / \ _ f ? +« p .( » +g 54 , . 

V(p,q-,p ,q) = c 3 e y ; e 1 J . (5.4) 


Oiy 2 

c 3 

e 2 

e 3 

1 

-0.1 

-0.920 

-3.225 

1 

-0.05 

-0.920 

-2.885 

1 

0.0 

-0.920 

-2.567 

1 

0.05 

-0.920 

-2.279 

1 

0.1 

-0.920 

-2.027 

2 

-0.1 

-0.474 

-2.571 

2 

-0.05 

-0.474 

-2.133 

2 

0.0 

-0.474 

-1.708 

2 

0.05 

-0.474 

-1.307 

2 

0.1 

-0.474 

-0.951 


Table 5.1: 3-body ground-state energies in fm -2 for different parameters chosen in 
Ref. [30]. a V 2 = 1 means that C\ = 12. fm -1 , c 2 = —12. fm -1 , o\ — 0.2 fm, and 
02 = 0.8 fm. a v2 = 2 means that c± = 0. fm -1 , C 2 = —2. fm -1 , and 02 = 0.8. fm. 


This potential has one very nice benefit of a clear separation of scales; each Gaussian 
has a well defined length-scale, a. With a choice of cr*, we can set up a “long-range” 
and “short-range” piece of the potential and carry out our J-shcll fits. The potential, 
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however is not ideal for our calculation for a number of reasons. The values of a, and 
Cj chosen in Ref. |3T7. 31'] require a very large number of points for accurate matrix 
elements when transforming back and forth between momentum representation and 
HO-basis, or SRG evolution (due to the diagonal never diminishing at high energies). 
Secondly, a Gaussian is not a good analogy for the long-ranged nature of the OPE 
potential. Lastly, as is shown in Ref P3], having a large ground state energy causes 
trouble if we need to evolve very low in A. 

Therefore, we will choose a better 2-body potential for studying universality in 
matrix elements. We still utilize a Gaussian for the short-ranged piece, but will use 
an exponential for the long-ranged piece. 


r 

V 2 (r) = -— e CT i + c 2 e _ T> 


CTly/n 


(5.5) 


We choose the constants such that there is some separation, r, above which the 
potential is almost entirely the long-ranged piece, and that the ground state binding 
energy is smaller than the A’s we wish to evolve to. For instance, for the values in 


Table C.l at r = 1.0, the contribution from the long-range exponential is —10.8 frn -2 , 
while the contribution from the short-range piece is only 1.9 x 10~ 9 fm -2 . Also, we 
must make sure that the short-range piece is non-negligible in the calculation, such 
that it acts as a hard core and is strongly repulsive at short-range. 

In Jacobi-momentum representation, the potential is: 

lV% &2 


~ , ci _ d(pV» 2 c 2 

V2{P,P) = - -?= e s + 


27ta/2 2tt °2 ip-p') 2 + 


(5.6) 


We then, analogous to realistic chiral potentials, multiply by a separable regulating 
function: 

r/w(*) a %-(£) 2 " 


V 2 (p,p')^V 2 (p,p')e-W e-W”. 


(5.7) 


75 








The following section will detail universality in the 2-body sector for our model 


potential, followed by our findings in the 3-body sector. 

5.2 2-body Sector 


For our study, we choose the values shown in Table C.l to construct our potential. 


Cl 

c 2 

O 1 

o\ 

A 

n 

48 fm -1 

-24 fm“ 2 

0.2 fm 

1.25 fm 

3.0 fm” 1 

4 


Table 5.2: Parameters for 2-body potential 


This potential has a single bound state with energy, Em = 0.402 fm -1 . The momen¬ 
tum regulator defines a well defined minimum for the UV cutoff from our H O-basis 
truncation. We choose b and iV max such that the truncation cutoff is above A = 4.0 
fm 1 , as the regulator at this value is negligible. We can choose a value for our 
coordinate-space truncation, L, such that the long-range part of the potential is neg¬ 
ligible. L — 10 fm is a safe estimate based on the Fourier transformation of the 
regulator. With these values in mind, we chose b = 2.0 fm and iV max = 120. The 
truncation cutoffs from this (roughly A ~ 7.8 fm -1 and L ~ 31.3 fm) are more than 
enough for accuracy in transforming between bases, allowing accuracy in evolution to 
lower values of SRG A. 


Fig. 5.1 shows from left to right: the model potential in momentum representation, 
in HO-basis, and the difference between the original and twice-transformed potential 
in momentum representation. We can see very clearly that above 4 fm -1 , momentum 
representation matrix elements of our potential are negligible, thus there should be 
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Figure 5.1: Left: Matrix elements V(p,p') of our model potential. Center: ffO-basis 
matrix elements of the same potential transformed into HO-basis. Right: Difference 
of matrix elements between the original potential and the twice-transformed potential, 
multiplied by 10'. 


no error in transforming our potential to HO basis. The cutoffs imposed at this 
momentum are e^f) 8 = 4.6 * 1CT 5 , thus we can safely say above 4 fm -1 , matrix 
elements are negligible with several digits of precision. Observing the HO-basis matrix 
elements, we can see that even an 7V max of 24 would be large enough to transform 
the potential between HO-basis and plane-wave bases. Finally, to explicitly show 
the difference in matrix elements of the original potential and the twice transformed 
potential in momentum-representation, the third picture shows these matrix elements 
multiplied by 10 7 such that they can be seen in the plot. This error in the 7th decimal 
is much more precise than we would expect from a typical realistic calculation, but 
must be kept low when looking for other, more substantial errors in calculation. We 
start with negligible errors in our basis-transformation and more points than necessary 
so that we can clearly identify errors that will arise from the SRG-evolution. 

The next step in our procedure is to generate a potential that will evolve to the 
same universal form in low-energy matrix elements as our model potential. Much like 
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Figure 5.2: (Left) for model, universal, and non-universal potential matrices in HO- 
basis; n = 2 i. (Middle) £(k) for model, universal, and non-universal potential matrices 
in momentum representation. (Right) Half-on-shell K matrix elements, K(k,k 0 ), for 
ko = 1.5 fm -1 . The value of £ is grid-dependent, thus the different scales for the first 
two plots. 


before, we construct this potential as the sum of the long-ranged exponential present 
in the model potential, but replace the short-ranged Gaussian with a sum of regulated 
5-shclls. Each shell has a different coupling which we fit to the lowest eigenvalues of the 
model potential. This potential conforms to the rules for universality we observed in 
Chapter [3] as it has the same explicit long-range potential and is low-energy eigenvalue 
equivalent to the model potential. We also create an eigenvalue-equivalent potential 
that does not have the explicit long-range piece of the potential (thus it will not evolve 
to a universal low-energy form) by fitting only h-shells without the long-range piece 
of the potential to the low-energy eigenvalues of the model potential. 


Figure R2 shows the eigenvalues in HO-basis and momentum representation and 
a vertical slice of the half-on-shell K matrix elements of the three potentials. Notice 
that the eigenvalues are very similar for all three potentials, but the low-energy HOS 
K-matrix elements for the non-universal potential are quite different from the model 
and universal potentials. At higher energies, the HOS K-matrix elements are all dif¬ 
ferent, as expected. The HOS K-matrix equivalence between the model and universal 
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Figure 5.3: (Left) HO-basis diagonal matrix elements of our three potentials, un¬ 
evolved (solid) and evolved in HO-basis to A = 2.0 fm -1 (dashed). (Right) Momen¬ 
tum representation diagonal matrix elements of our three potentials, unevolved (solid) 
evolved in momentum representation to A = 2.0 fm -1 (dashed). 


potentials seems to end approximately above 1.5 fm -1 . While we are still concerned 
with 2-body potentials, we have the luxury of being able to plot momentum represen¬ 
tation eigenvalues and HOS K matrix elements, so we have done so to double check 
that the procedure is the same as before. We can see that the model and universal 
potentials exhibit the criteria to evolve to a universal form, and the nonuniversal, yet 
low-energy eigenvalue-equivalent potential does not. 

Now, we can evolve the potential matrix elements in momentum-representation 
and HO-basis. Observing Figs. 5.3|5.6 , we see the emergence of universal matrix 
elements. These figures show potential matrix elements of the three potentials in 
(left) HO-basis and (right) momentum-representation both unevolved (solid lines) 
and evolved (dashed lines). The potential designed to be non-universal, as planned, 
does not exhibit a universal low-energy form. We can see this in Figs. 5.3||5.4 Notice 
that below 1.5 fm -1 , the momentum-representation matrix elements evolved to A 
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Figure 5.4: (Left) HO-basis off-diagonal matrix elements of our three potentials, un¬ 
evolved (solid) and evolved in ffO-basis to A = 2.0 fm -1 (dashed). (Right) Momentum 
representation off-diagonal matrix elements of our three potentials, unevolved (solid) 
evolved in momentum representation to A = 2.0 fm -1 (dashed). 



Figure 5.5: (Left) HO-basis diagonal matrix elements of our three potentials, un¬ 
evolved (solid) and evolved in HO-basis to A = 1.0 fm -1 (dashed). (Right) Momen¬ 
tum representation diagonal matrix elements of our three potentials, unevolved (solid) 
evolved in momentum representation to A = 1.0 fm -1 (dashed). 
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Figure 5.6: (Left) HO-basis off-diagonal matrix elements of our three potentials, un¬ 
evolved (solid) and evolved in ffO-basis to A = 1.0 fm -1 (dashed). (Right) Momentum 
representation off-diagonal matrix elements of our three potentials, unevolved (solid) 
evolved in momentum representation to A = 1.0 fm -1 (dashed). 


= 2.0 fm -1 of the model and universal potentials (black-dashed and blue-dashed 
lines respectively) fall on top of each other, while the non-universal potential matrix 
elements (red-dashed) are very different. Viewing the ffO-basis matrix elements at 
this A is inconclusive, as only the first two matrix elements in ffO-basis are the same 
between the universal potentials, while the rest differ. 

Evolving further to A = 1.0 fm -1 (Figs. 5.5|5.6 ), we see that more matrix elements 
align in ffO-basis, but also that the liigh-n diagonal matrix elements are no longer 
zero. This is also true for the (not-shown) near-diagonal elements. This means that 
transformation between bases will not be as accurate as with the unevolved potentials. 
As A decreases further, oscillator errors will appear in the high-momentum matrix 
elements of the potentials in momentum representation. 

We can see that momentum representation is ideal for viewing matrix elements and 
searching for a low-momentum universal form. A momentum scale appears around 
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1.5 fm -1 , below which matrix elements are universal, and it is easier to identify this 
universal form at A > 1.0 fm -1 in momentum representation than in HO-basis. It 
may be possible to identify universality in HO-basis matrix elements, but as we evolve 
to lower A, we accumulate SRG errors in HO-basis and potentially run closer to the 
binding momentum which creates pathologies in the evolution with ^rel- In the 3- 
body problem, we cannot evolve in momentum representation, nor can we evolve 
to as low a A. Our strategy for the 3-body problem will be to evolve in HO-basis 
and observe matrix elements in HO-basis and then those transformed into momentum 
representation. We have chosen our HO-basis parameters to have much more precision 
than required in the 2-body sector, such that we can isolate imprecisions that appear 
in the 3-body problem. 

5.3 3-body HO-basis evolution 

After embedding the 2-body potentials of the previous section, we can calculate 
bound state energies. The model 2-body potential has a 3-body ground state energy of 
0.991 fm 1 , the universal potential a ground state energy of 0.989 fm -1 , and the non- 
universal potential 0.961 fm -1 . Each potential has a number of three-body excited 
bound states. None of these ground state energies need be identical; as has been shown 
in Ref. [55]; different universal 2-body potentials may yield different 3-body ground 
state energies. Realistic 2-body potentials do not predict exactly the same triton 
binding energy, and the predictions do not reproduce the experimental value without 
a 3-body force [6]. The predictions should not be off by an order-of-magnitude, 
however, as this could be evidence of strong many-body forces that break the many- 
body hierarchy. 


82 


When evolving in HO-basis, it is standard (for both realistic and other model 
potentials) to build the entire 3-body Hamiltonian and then plug it into the SRG 
flow equation. The method of Refs. [ST, T6] [TIJ to subtract out the disconnected 
diagrams in the 3-body evolution is naturally defined in a 3-body plane-wave basis, 
but much more difficult, if possible, in HO-basis. After evolution, we can isolate the 
3-body piece by evolving the 2-body potential in a 2-body basis, embedding it into a 
3-body basis, and then subtracting it from the total 3-body potential (embedded 2- 
body plus explicit 3-body) evolved in the 3-body basis. This evolution is unitary and 
will preserve the binding energies at any A up to numerical errors in the differential 
equation solver. 

Once the evolved 3-body potential (which can contain both explicit and induced 
pieces) is isolated, transforming back to momentum representation is straightforward. 
When viewing potential matrix elements of the evolved 3-body potential in momen¬ 
tum representation, oscillations indicative of basis transformation errors appear at A 
values much higher than expected. In fact, the oscillations appear even at A values 
so high that no discernible changes occur in the 2-body evolution (or in an explicit 
3-body potential we could include). Our method accurately transforms unevolved 


3-body potentials (Eq. (5.4)) from HO-basis to momentum representation, thus it is 
clear that the evolution causes the oscillations in momentum representation of evolved 


3-body potentials (see Fig. 5.7). 


Fig. |5.8| shows momentum representation matrix elements for the induced 3-body 
potential. We select matrix elements with zero spectator momentum and see oscilla¬ 
tory behavior in the low pair-momentum elements which increases in magnitude as 
A is decreased. This oscillation is similar to what we observe for a basis truncation 
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Figure 5.7: Momentum representation matrix elements diagonal in pair momentum 
with q=q’=0 of the explicit 3-body potential, V(p,p’;q,q’). The dashed red line 
is in momentum representation where the potential is initially defined. The solid 
blue line is the potential transformed into HO-basis and combined with the 2-body 
embedded potential, which is then removed and then the matrix is transformed back 
to momentum representation. The maximum difference in any matrix element is 
7.61xlCT 8 . The 3-body potential used is from Eq. (5.4) with C 3 = 1.0. 


84 













Figure 5.8: Momentum representation matrix elements diagonal in pair momentum 
with q=q’=0 of the induced 3-body potential, V(p,p’;q,q’), (left) evolved to A = 14 
fm _1 and (right) evolved to A = 1.5 fm _1 . Notice the magnitude of the oscillation in 
these low-momentum matrix elements increases as A decreases. The oscillations ap¬ 
pear due to truncation of nonzero high-n HO-basis matrix elements, and thus appear 
as high-frequency oscillations. 


error in the 2-body sector; appearing as high-frequency oscillations due to omitted 
high-n oscillator wave-functions. Furthermore, the oscillatory behavior appears at 
much higher A than we expect, based on the evolution in the 2-body sector (or lack 
thereof). For this reason, the oscillatory behavior most likely is not representative of 
the induced 3-body potential, but rather a truncation error on the matrix elements 
of the full, evolved 3-body potential. 

For further evidence that the momentum-representation oscillatory matrix ele¬ 
ments of the induced 3-body potential are a sign of basis-truncation errors in the 
matrix elements, we can examine the HO-basis matrix elements of the embedded, 


unevolved 2-body potential. Fig. |5.9| shows a selection of off-diagonal potential ma¬ 
trix elements of the unevolved embedded 2-body potential. We can see that the 
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Figure 5.9: HO-basis matrix elements of the unevolved embedded 2-body potential 
as a function of n p with n' p = 0,n q = 0 and (left) n' q = 60, (middle) n' q = 90, (right) 
n' q = 120. Notice the shape and decreasing magnitude of each sample of matrix 
elements as n' q increases (the vertical scale is different in each plot to better resolve 
shapes). 


matrix elements for n' = 60 and 90 are roughly localized just above the point where 
n p + n q = n' p + n' q . Therefore, when we look at the extreme off-diagonal matrix el¬ 
ements (n' q = 120), we can see that they are increasing up to the basis truncation 
(iV max = 120), and using the shape when n' q = 60 and 90 as a guide, we presume 
that a higher truncation would produce nonzero matrix elements at n p values slightly 
greater than 120. Much like the oscillatory truncation errors that appear in the 2- 
body sector, these truncated 3-body matrix elements produce oscillatory errors when 
evolving in the 3-body sector. Unlike the 2-body errors, these truncated matrix ele¬ 
ments are unavoidable in the total n truncation scheme we use. 

A truncation on each individual n, specifically n p much greater than n q , might 
solve the problem of truncating non-negligible matrix elements from the embedded 
2 -body potential, but the total n truncation is required for a unitary symmetrizer. 
If the symmetrizer is not unitary, there will be errors in the symmetrized basis. 
Furthermore, a much more inhibiting problem with a non-unitary symmetrizer is that 
the eigenvalues will no longer be either zero or one, but somewhere in between. This 

86 









means that at best one must choose only nearly symmetric states (eigenvalue close to 
1 ), and throw out barely symmetric states (eigenvalue close to 0), and at worst states 
may have eigenvalues in between such that the process cannot differentiate symmetric 
states from the rest at all. Lastly, some of the matrix elements we would omit in the 
symmetrizer with a truncation on each n would be the very same far off-diagonal 
matrix elements that are causing the error in SRG evolution, thus the scheme with 
the least problems is a total-n truncation. 

At the heart of the problem is the 2-body potential between one of the pair 
particles and the spectator. The momentum entering is a combination of the pair 
and spectator momentum, and thus because the spectator momentum (or oscillator 
number) can be anything, the embedded potential will have coupling between low 
and high pair momentum (or oscillator number) matrix elements. In the unevolved 
potential, these far off-diagonal matrix elements play a small role, thus we can find ac¬ 
curate binding energies, but as the potential evolves in HO-basis, the errors in matrix 
elements accumulate with each step in the ordinary-differential-equation solver. 

To further clarify these errors as an error in the HO-basis SRG evolution and not 
just as a part of a pathological induced 3-body potential, we can examine the 3-body 
evolution equation only using the first term, quadratic in T. This term does not 
induce any 3-body potential, and without an explicit 3-body potential it is identical 
to the first term in the 2-body evolution equations for each 2-body potential. We can 
see each 2-body evolution equation with just the leading-order term looks like: 
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Figure 5.10: Error in momentum representation matrix elements of the embedded 
2-body potential evolved with only the quadratic in T re i term (q — q' — 0 and p = p'). 
(left) evolved to A = 14 fm _1 and (right) evolved to A = 1.5 fm -1 . Notice the 
magnitude of the oscillation in these low-momentum matrix elements increases as A 
decreases. 


The 3-body evolution equation with just the leading-order term, 


dVu 


EVT = E[[ T rel’ V ii], T rel] • 

i¥=j *7 


can be simplified using: 


[TrehVij] — [Tij, Vij] 

[ [Tij ,Vij],T re i] = [ [ Tij , Vi j ], Tij ], 


(5.9) 


(5.10) 

(5.11) 


which are easily worked out in either HO-basis or momentum representation. We 
then see that the 3-body leading-order evolution equation is simply a sum of each 
of the three 2-body terms. Therefore, embedding the 2-body potential into the 
energy-truncated symmetric basis and evolving with the 3-body leading-order equa¬ 
tion should be identical to evolving the 2-body potential with the leading-order term 












in the 2-body basis, then embedding into the energy-truncated symmetric basis. Ob¬ 


serving Fig. 5.10 we can clearly see the same oscillatory errors appearing in the mo¬ 
mentum representation matrix elements which get worse as A decreases. Therefore 
it is clear that there is an error in the HO-basis 3-body evolution that increases as A 
decreases. 

It is important to note that this is an error in the matrix elements of the 3-body 
evolved potential. It will not appear as an error in the eigenvalues or any 3-body 
observables calculated in HO-basis, as the SRG is a unitary transformation. If the HO- 
basis potential was transformed to momentum representation, and a calculation done 
in momentum representation, errors in calculated observables would likely appear. 
As adding additional bodies increasingly adds importance to off-diagonal potential 
matrix elements through additional momentum integrations {56], we expect to see 
the errors in potential matrix elements translate to many-body errors in observables, 
increasing with particle number. Most apparent, however is that our prescription for 
observing momentum representation matrix elements in a search for universal form 
after evolution is impossible with the oscillatory errors arising from the evolution. 

The errors in HO-basis evolution do, however, decrease with a larger basis size. 


We see in Fig. 5.9 that the far-off diagonal matrix elements decrease in magnitude 
as total N increases. For an accurate HO-basis evolution, we see that A max must be 
taken to be much larger than is necessary to reproduce 3-body observables, and for a 
desired accuracy at a given A, it may be possible to simply make N max large enough. 
It may also be possible to bypass the problem using the evolution method of EH in 
a plane-wave basis, where it is well defined. It may also be possible to subtract out 
the disconnected diagrams in HO-basis and evolve in an analogous way. 




5.4 3-body Recap 


We began this chapter confirming ground-state energies for a previously used 1- 
D potential for spinless bosons [301 Eg. We introduced a 2-body potential more 
beneficial to our study because its long-range interaction is not as confined as a 
Gaussian, it requires smaller basis size for accurate calculations, and it is easier to 
fit with our 5-shcll procedure. We confirmed that this potential follows the same 
rules for universality in evolved potential matrix elements as the realistic potentials 
of Chapter [3j This potential evolves accurately to A = 2.0 fm -1 , but continuing the 
evolution to A = 1.0 fm -1 shows signs of truncation errors. 

Continuing into the 3-body sector, we see that oscillatory evolution errors appear 
at much higher A than in the 2-body sector. Small errors appear at A = 14 fm -1 , and 
evolving further to A = 1.5 fm -1 increases the errors to be of comparable magnitude 
to the unevolved potential. We identify the source of these errors as far-off diagonal 
matrix elements arising from the spectator particle. To be sure that these oscillations 
are SRG evolution errors, we show the difference in matrix elements evolved with 


the term quadratic in T of Eq. (5.8), which should be the same if calculated in the 
2-body sector then embedded into 3-body or calculated in the 3-body sector. Fixing 
this SRG error in 3-body SRG evolution is a critical step in more precise Hamiltonians 
for many-body calculations. 
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Chapter 6: CONCLUSION 


6.1 Recap 

The modern precision 2- and 3-nucleon potentials are the product of decades of 
discovery and refinement, but their success in few-body systems is not yet mirrored 
in many-body calculations. Even the soft yEFT potentials possess too strong a cou¬ 
pling between relevant low-energy and high-energy potential matrix elements, thus 
decoupling these matrix elements through SRG evolution is a critical step for solving 
many-body problems. In an effort to improve 2-body potentials, we can use univer¬ 
sality in evolved potential matrix elements as a guideline to simplify required aspects 
of the potential. We have seen in chapter [3] that for matrix elements of different 
2 -body potentials to evolve to the same universal form in a low-energy region, their 
unevolved HOS T-matrix must be equivalent in the same energy region, and they 
must have the same binding energies (if any). Potentials possessing the same explicit 
long-range terms exhibit the HOS T-matrix equivalence required for universal low- 
energy potential matrix elements. One may also utilize this in reverse; understanding 
the effective range of an operator based on the effect in the HOS T-matrix elements. 

The necessity of an explicit 3-body potential for accurate 3- and many-body calcu¬ 
lations means that it is important to extend our understanding of universality into the 
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3-body sector. The complexity of adding an explicit 3-body interaction that comple¬ 
ments the 2-body potential has limited our methods for solving few- and many-body 
nuclear problems. For this reason, evolution of potentials with explicit and induced 
3-body forces for use in many-body calculations has been restricted to HO-basis with 
limited maximum oscillator number. There is far less understanding of evolution 
precision in 3-body spaces than in 2-body, and explicitly observing differences in 
momentum and HO-evolution is important. 

We have observed in chapter [5] that truncation errors imposed by the HO-basis 
are manageable in the 2-body sector; one has the ability to set momentum cutoffs and 
mesh sizes sufficiently high to achieve desired precision. In 2-body potential matrix 
elements, it is possible to only truncate negligible matrix elements. We observe that 
in the 3-body sector, the spectator particle adds strength to far off-diagonal matrix 
elements in the embedded 2-body potential. The strength of these far off-diagonal 
elements decreases with our truncation of oscillator number, but at too slow a rate, 
which inhibits computation. Errors in the 3-body observables from truncation of these 
non-zero matrix elements in the 3-body sector can be overcome with sufficiently large 
mesh size, however due to the nature of HO-basis SRG evolution with the T generator, 
the errors compound and grow as A decreases. We observe this as oscillatory errors 
in momentum-representation and differences in the quadratic-in-T term evolution of 
2-body potential matrix elements in the 2-body space and 3-body space. 

More precise evolution of embedded 2-body potentials in HO-basis is required 
before universality in evolved potential matrix elements in the 3-body sector can be 
studied in a way analogous to chapter [3] Furthermore, these errors in potential matrix 
elements arising from HO-basis SRG evolution do not affect the 3-body observables, 
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but they could be a source of error in many-body calculations. As an increasingly 
important focus of low-energy nuclear physics is estimating theoretical error bars, 
a better understanding of the HO-basis evolution errors is important, and we have 
identified an apparent source of this evolution error. 

6.2 Looking Forward 

The natural next step in this study will benefit greatly from 3-body evolution in 
momentum representation. Once precision evolution in a momentum basis (plane- 
wave Jacobi pj6] or otherwise Ed) is accomplished, it will be straight-forward to 
compare with HO-evolution and identify the extent of the differences. Some many- 
body methods will still require the computational efficiency of HO-basis, but if the 
errors of HO-basis evolution are indeed large, a momentum-basis evolution followed 
by transformation to HO-basis would create better input Hamiltonians for many- 
body calculations. Furthermore, we saw that observing universality in evolved matrix 
elements is much more straight-forward in momentum representation, thus a study 
of 3-body universality would be much clearer if precise momentum representation 
evolution was possible. Should 4-body forces be recognized as important contributions 
to the input Hamiltonians, HO-basis evolution will certainly be much simpler to 
achieve than momentum representation evolution, and the lessons learned about the 
differences in the two evolution methods in the 3-body sector may provide valuable 
insight into the errors in 4-body evolution, where direct examination of the differences 
may not yet be possible. 

Because of the difficulty of generating explicit 3-body potentials, far fewer of them 
exist than 2-body potentials. In principle, each 2-body potential may have an infinite 
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set of complimentary 3-body potentials which together accurately reproduce observ¬ 
ables in the 3-body sector. Understanding universality in the 3-body potentials may 
serve as an important simplification to generating them and may identify features 
required to improve accuracy of few-body input Hamiltonians in many-body calcula¬ 
tions. 
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Appendix A: ISSP DERIVATION 


The following is the derivation of the ISSP from Ref. ra- Starting from the 
Lippmann-Schwinger equation (LSE) in partial-wave momentum representation (omit¬ 
ting Z-indecies), 


(ik\K(E)\k ') = {k\V\k') + -P 

7T 


dpp2 (k\V\p)(p\K(E)\k') 


(A.l) 


we can input a separable potential, 


(k\V\k') = av(k)is(k'). 


This simplifies the LSE to, 


(k\K(E)\k f ) = 


cr v(k)v(k') 
D(y/E) ’ 


D(x) = 1 - -P f dp p 2 V ^ 

7T J X z — p z 


We can define, 


D + (x) = lim (D(x + ie)). 

K y e —>0 V 


Evaluating this limit, we find, 


D + (x) = D(x ) + iov 2 (x)x. 


Combining this with equations A.4 and 3.8[ we see, 

Im[D + (x)\ 


tan(— S(k)) = 


Re[D+(x )] ' 


(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 
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It follows to write, 


D + (x) = \D + (x)\e~ i5{k \ 

(A.8) 

ln(D + (x)) = ln\D + (x)\ —i8(k). 

(A.9) 


This final equations shows that the imaginary part of ln(D + (x )) is determined by the 
phase shifts. 

The next step requires making a contour integral. Given a function, f(z) analytic 
in a closed contour, C, containing point, z, we shall use Cauchy’s theorem to write, 

f(z) = — & dx^-. (A.10) 

2tt i J x — z 

If we choose the contour including the real axis and a semicircle in the upper-half 
plane (UHP), and the function vanishes as \z\ —> oo in the UHP, then adding ie to z 
to guarantee we remain in the UHP this becomes, 


f(z + ie) 


M . 

2m J_ 00 x — z — ie 


(A.11) 


Taking the real part, we see 


1 

Re[f(z + ie )] = - 

TT 


dx 


Im[f(x) 


x — z 


(A.12) 


We know that the potentials we consider fall off as k —> oo, thus we can omit 
the semicircle part of the contour. We also know that a bound state will produce a 
simple pole in the K-matrix, which is caused by a zero in D + (\/E). We can modify 
D + (x ) in a way that will remove the zeros in the UHP for our contour integration, 
and also preserve the asymptotic behaviour of the function. Introducing Kb = \/E B , 
this modification is 

xv h -A- i K r> 

(A.13) 


D+(k) = P~^S-D + (k). 


k — ik b 
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Now we can relate the real part of ln[D + {k)\ to an integral over the imaginary part, 
which is simply the phase-shifts. We start by noticing that 


Re[ln(D + (k))] = hi(D + (k )), 
Im[ln(D + (k))] = — iS(k ) + ln{ 


k + m B 

k — kg' 


(A.14) 
(A.15) 


Now, using A.12, we see 


ln\D + (k))\ = -A (k) + —P 

7rr 

\ / j \ i „ r dk 

= * p L-* 


°° dk' i ,k + inBs 

iri J_ OQ k'~k v k — ins h 
dk'5(k ) 
k' 


(A.16) 
(A.17) 


The second term (integral) in A. 16 is evaluated as ln( k ^ B ). It follows, then that 

(A.18) 


7 2 2 

= k + K B „-A(k)-iS(k) 


k 2 


Now, we simply equate the imaginary part of this equation, and equation |A.6[ to find 


au (k) = 


k 2 + k% sin (k) A(k) 


k 2 


(A.19) 
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Appendix B: FURTHER HARMONIC OSCILLATOR 

BASIS DERIVATION 


In this appendix we show the derivation of T re i in HO-basis for 2- and 3-bodies 
and the derivation of the transition bracket. 

To calculate T^ 2 \ the 2-body relative kinetic energy operator, we utilize ladder 
operators m ■ 


(n\T^\m) = 

p2 

(n\ —m) 
x 1 2 m 1 ' 

(B.l) 

P = 


(B.2) 

(n T ( 2 ^ \m) = 

x 1 2 m 17 

(B.3) 

(n T 1 - 2 ^ m) = 

(n-— (aW — a*a — <W + aa) m) 

' 4b 2 

(B.4) 


(n|T ( 2 )|m) = — — [ yjm + 1 yjm + 2(n\m + 2) 

— ( y/myfm + y/m + 1 y/m + l)(n|m) 

+a Jn + 1 y/n + 2 (n + 2\m)\ (B.5) 

{n\T^\m) = ~^[ y/m + ly/m + 26 n>m+2 

(2n T l)5n,m 

+a/?t. + 1 y/n + 2^ n ._|_2 i m]. (B.6) 
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To calculate T (3 ^ we take advantage of the fact that | n p ) and \n q ) are in orthogonal 
spaces, thus: 


(n p n q \T^ \m p m 

\ / A W i \ 

(B.7) 

(n p n q \T^ \m p m 

P 2 Q 2 

q ) = (n p \(n q \— + — \m p )\m q ) 

(B.8) 

'n p n q \T (2) \m p m q ) = 

P 2 Q 2 

{n p \ — \m p ){n q \m q ) + (n q \ — \m q )(n p \m p ) 

(B.9) 

[n p n q T (2) \m p m q ) = 

(n p T (2) \m p )5 nqtTnq + (n q \T {2] \m q )d np:mp . 

(B.10) 


Lastly, we must calculate matrix elements of the transition bracket, {n p n q \P 2 z\n' p n'^) 
(n p n q \n' p n' q ) 3 . The Jacobi momentum with particles 2 and 3 exchanged is readily cal¬ 
culated from Jacobi momentum definitions. After some algebra, 


P = 


Q = 


1 , V3 

2 P + T 9 ’ 

V3 1 

irP-nQ- 


(B.ll) 

(B.12) 


We can define ladder operators for each oscillator basis state, r/,-" 1 10) = -?=\ rii), which 
transform exactly the same as momenta: 


% 


i , V3 
2 T,P + ~2 
a/3 


~Vp 


1 

2 


Vq, 

(B.13) 

Vq- 

(B.14) 
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Using these definitions, 


{n'n'\n p n q ) 3 = 


(otfSN^lo) 


y/ri p lri q lnplnq 


n’ p \n' q \n p \n q \ 2 


lrJ / 

T 




/„*, , T <°l E ) [k] (Bt “ fc) [^d 

vHKSS 1 f^\kj 2 2 


9 / / 

n. 


x 


Emi^bhrf 

j=o v J ' 

X^np^lQ) 


EE 


^In^npln,! “ “ V ^ / V i 


K 


X 77/p! Tlq \8 n / k-\-n' —j,rip &k-\- j,n q 


in n in„i sr^ fn f p 


Lp.lLq. \ a 

7 / I 77 / ! / 


n' 


1 ) r ^ I (Wp+Wg—2fe) r (w^—n 0 +2fc) 

n' p vri' q \ ^ \ k ) \n n - k) 2 2 


x6(n q — n q + k)6(n q — k ) 

X ( —l)«9-^ 


(B.15) 


(B.16) 


(B.17) 


(B.18) 


(B.19) 


The third line utilizes the binomial theorem and the fourth line balances the ladder 
operators. 

From this we can see why we need a truncation in total n p + n g . If the truncation 
is made individually on each, we can examine the lowering operators in the previous 
equations. They enter in powers, [} 2 q p + ^ Vq] n ' p [ 2 Vp ~ \ r lq\ n ' q ■ We see that each 
lowering operator can enter n' p + n' q times, connecting our originally truncated states 
to states with oscillator number above our truncation. This makes our symmetrizer 
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non-unitary, which causes eigenvalues to no longer be distinctly 1 or 0. Rearranging 
the states and truncating on total oscillator number, however, solves this problem. 
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Appendix C: TABLE OF ABBREVIATIONS 


AVI 8 

Argonne Ui 8 

HO 

harmonic oscillator 

HOS 

half on shell 

IR 

infrared 

ISSP 

inverse scattering separable potential 

LS 

Lippmann- Schwinger 

OPE 

one-pion exchange 

QCD 

quantum chromodynamics 

SRG 

similarity renormalization group 

uv 

ultraviolet 

xEFT 

chiral effective field theory 


Table C.l: Table of abbreviations 
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